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PROBABILITY DISTRIBUTIONS ASSOCIATED WITH 
DISTINCT HITS ON TARGETS 


D. A. SPROTT 
DEPARTMENT OF PSYCHIATRY 
UNIVERSITY OF TORONTO, TORONTO, CANADA 


Some probability distributions connected with distinct hits on targets, 
using two different firing schemes, are developed. It is assumed that any 
shot has a probability p, not necessarily unity, of hitting the target at 
which it was aimed. The development uses a well-known expression for 
the probability that exactly ¢ of N possible events occur simultaneously. 
Some of the formulae developed here include as special cases the proba- 
bilities derived separately and by more complicated arguments in papers 
by N. Rashevsky (Bull. Math. Biophysics, 17, 45-50, 1955) and A. Rapo- 
port (Bull. Math. Biophysics, 13, 183-38, 1951). 


The probability of having exactly m empty cells when s objects 
are distributed at random to N cells 


N-m 


Prcm| s)=() ¥* cay (“5 ")a- im 
j=0 


is well-known and has proved to be very useful. Recently, modifi- 
cations of this formula have been used in certain biological problems 
(Rapoport, 1951; Rashevsky, 1955). Rapoport considers the proba- 
bility of hitting & new targets with s shots, having hit m targets 
with s shots previously, where there are N targets and it is assumed 
that each shot hits some target. Rashevsky considers the problem 
of assigning biological functions to cells. In this note are con- 
sidered various probability distributions when there is a proba- 
bility p #1 of a shot hitting its target. Such distributions may also 
arise when we sample populations, divided by some criterion into 
groups, for some characteristic which occurs with a frequency p. 
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Probability distribution of the number of targets hit. Suppose s 
shots are fired randomly at N targets and that there is a probability 
p that any specified shot hits its target. Let 7 specified targets 
have s shots fired randomly at them. The probability of all s shots 
missing the j targets is (1 - jp/N)*. Thus, in the notation of Fel- 
ler (1950, p. 61), 


$= 0101 - ip/N) = (") (1 = jp/N)*. 


Hence the probability that exactly ¢ targets are missed is the prob- 
ability that exactly N - ¢ targets are hit, which is 


Pr(W = tls) = 7 (2) 57) 84 


=O 


N-t 


=> (-1)/ ("5 i? ) {1 —(j + é)p/N)§ 


j=0 


ee iS ay ("5 ‘) (1 - (i+ Op/Mi. 


The probability that exactly m targets are hit is 


ErGn ley (*) 5. (1)! (7) [1 (N= 4 f\peNien 


29 


If p = 1, this gives the usual formula 


Pr(m|s) = ie De 2B (-1)/ (j (m— j)*/N*. 


Various forms of this formula occur in DeMoivre (1718, pp. 109- 
110, problem 39), Whitworth (1934, pp. 87-88), Stevens (1937), 
Fisher (1949, pp. 11-12), Feller (loc. cit., p. 69); it is also con- 
sidered by Euler, Laplace, Mallet, and Trembly (see Todhunter, 
1865, pp. 159, 250-257). 

Rashevsky (loc. cit.) defines RN to be the number of different 
ways that s objects (biological functions) can be distributed to N 
non-empty classes. Since, if the classes were boxes or cells they 
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could be permuted in N! ways, and since there are N° ways of 
distributing s objects to N cells without restrictions, we have the 
obvious relation 


N 
NIRN/NS = Pr(m=N|s, p=1) = * (-1)! (7) (1 - j/N)*. 
AD t 
Hence 
* N 
RN = -1)/ (". — )*/N! 
B= ay (Y) ww an 
j=0 
as given by Rashevsky (loc. cit.). 
Suppose that wu additional shots are fired; as before, the proba- 
bility that j specified targets are missed is (1- jp/N)”. Thus, 
since only the N — m targets missed previously are being considered, 


8, = ("5") a ipmy. 


The probability that exactly vw of these N — m targets are missed, 
which is the probability that exactly N — m - v of them are hit, is 


N-m-v 
Pr(N- m-—2|v) = oH De cay (MPN (ws Dern. 


v 
j=0 


The probability that exactly & additional targets are hit is 


k 
Pr(k|u) = 3") ye 1)! () [1-(N-m-—k+j)p/N]*. 


j=0 


The result for u= 8s, p = 1, is given by Rapoport (loc. cit.) 


N-m\ < ; (k 
Pr(tlu= 9, p= 1) ( k ae (j) me NONE. 


j=0 
To find the mean and variance of k, let 


X, = 1 if the ¢th target is hit, 
= 0 otherwise. 


Then 
E(X,) = Pr(X, = 1)=1-(1-p/N)’. 
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The total number of targets hit is 


thus the expected number of targets hit is 


N-m 
E(Ty__, | = ye E(X,) =(N-m)[1-(1-p/N)"). 


im 


The variance of the number of targets hit can be found in a similar 
way. 


E(X,X,) = Pr(X,X, = 1) =1-2(1 - p/N)* + (1-2 /N)". 
Thus 
var (Ty |u) = var ). X, =(N- m) var X, + 


(N — m)(N — m - 1) cov (X;, X,) 
= (N -— m) E(X?) -— (N - m) E?(X,) + 
(N-m)(N-m- 1) E(X,X,) - 
(N — m)(N ~ m-1)E?(X,) 
= (N — m)(1 — p/N)* — (N — m)?(1 — p/N)?2* 4+ 
(N —m)(N-m-1)(1-2p/N)*. 
The corresponding formulae for Pr(m | s) are obviously 
N{1 -(1- p/N)*], 
and 
N(1 — p/N)§ — N7(1 — p/N)?8 + N(N-1)(1-2p/N). 
Probability distribution of the number of shots required to hit 


exactly m targets, Pr(s|m). The probability that exactly m — 1 tar- 
gets are hit with s - 1 shots is 


N m-1 re 
Fs oe cay ( j )(1-@W-mere p/m, 


f=0 


The probability that one more target is hit with one more shot is 
(N-m-+1)p/N. Therefore the probability that exactly s shots are 
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required to hit exactly m targets is 


Pr(s|m) = p ee) i (-1)) (" F ‘ [1-(N-m4+1+j)p/NIS"?. 


j=0 


The probability that exactly s shots are required to hit all targetsis 


= Neeet 
Pr(s|m=N) =p ay ay ( ; ) t= Ge ayp/nye 


Oo 


The mean and variance of the number of shots required to hit 
exactly m targets can be found by the use of random variables. Let 
X, be the number of shots required to hit the next different target 
after hitting & different targets. 


Pr(X, = jf) =(1-(N-4)p/N)"! (N- 4) D/N, 


since the first 7 - 1 shots must hit targets that have been hit pre- 
viously and the last shot must hit a new target. Then 


m—1 
= 3 X, 
k=0 


is the total number of shots required to hit exactly m targets. 


E(X,) =), j{1-(N- ®)p/NII~ (N - k)p/N = N/p(N - &). 
jai 
Similarly, 
var X, = E(X?) - E?(X,) = N2/p? (N — k)? — N/p(N — k). 


Hence, 
m-1 
E(s|m) = )\E(X,) =(N/p) )| /(N- 4%) 
k =0 


and 


=o 


m-1 
var (s|m) = (N/p)? a; [1/(N — k)]? - (N/p) iE 1/(N - k). 


m 
k=0 k=0 
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In particular, 


N 
E(s|m=N) =(N/p) e 1/j —» (N/p)(y + log, N) where y = 
je=1 
.577216 (Euler’s constant), 
and 


N N 
var (s|m = N) =(N/p)? Dis (1/7)? - (N/p) ». iff 
je=1 i= 
(N 7)7/6 p? — (N/p)(y + log, N). 
Even distribution. Assume that s shots are fired at N targets in 
such a way that no shot is fired at a target that has been hit pre- 
viously. Thus, exactly m targets will be hit if exactly m of the s 
shots are hits; the probability of this is 


Pr(m|s) = () p™ q®™ where (¢ =1-p). 


m 


The mean of m is sp, and the variance is spq. 
Proceeding as before, the distribution Pr(s | m) can be found; the 
probability that exactly m - 1 targets are hit with s —- 1 shots is 


a= 1 m-1.s—m 
Pee fs q ; 


The probability that another target is hit by another shot is p, and 
so the probability that exactly m targets are hit with exactly s 
shots is 


: -1 m ~s—m 
Pr(o|m)= (2-4) p ee 


The generating function is 


f(t) = 2: f P ) O g*T eee Be Oa ed 


s=m 
Therefore the mean of s is 


E(s|m) = f’ (1) = m/p; 


N 
E(s|m = N) = N/p<(N/P) % (1/7), 


al 
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which was the expected number of shots required to hit all targets 
when the shots were randomly fired. The variance of s is 


var (s|m) = f’*(1) + f*(1) - f’(1)? = mq/p?. 


Appendiz. The previous expressions for the probability of exactly 
m hits are obtained from the general formula for the probability that 
exactly ¢ of N possible events occur simultaneously, 


N-t 5 
Pr(t) = oa (-1) ee Srey (1) 


j=0 


where S, is the sum of the (?) k th-order probabilities; a typical 


kth-order probability is p(z,,..., z,), which is the probability of 
the simultaneous occurrence of the events ty» tgye++,%,- In the 
application to hits on targets, ¢, is the event that target i, is 
missed by all s shots. 

A direct proof of the above formula for Pr(¢) is given by Feller 
(loc. cit., pp. 64-65), and a direct proof of a similar expression is 
given by Whitworth (loc. cit., pp. 68-71). It is simpler, however, 
to proceed indirectly by applying the following two theorems. 


Theorem 1. 


Xt 
s,-)° (*) pre, (2) 


k=j 


where the symbols are defined above. 

Proof. S, is the sum of probabilities of the form p(Z,, %,+++5 
ts woe, thy U4 40ee*s ty), Where 2, denotes the non-occurrence of 
the eventi,, and&>j. The probability p(7,, t,5+++5 jreees Upp 
Coeur. <p ba) occurs in the expression for Pr(s) if and only if 
s =k, and then occurs with coefficient unity. Therefore it occurs 


on the R.H.S. of (2) with coefficient (i) Since S, is the sum of 


the probabilities of the occurrence of j events, the coefficient of 
P(E qr Saye ees Biren Uys there? ty) in S, is the number of sub- 
sets of j events that can be formed from the & events tis tgyeees Ups 


that is, (i) . Since & can range from j to N, the theorem is proved. 
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Theorem 2. (The inversion formula.) 


If 
N 
fe — 
nS O29 
kaj 
then 
= j-c jm 
Wey (23) cay ii 


juk 
Proof. Substitute the expression for X, into the expression for Y, 
and invert the order of summation, noting that () =0 when a< b 
and that therefore 


x 3 ok 
ig : (2 ‘=o 


unless 7 = &, in which case the sum is unity. 
To prove (1), apply Theorem 2 to Theorem 1, with c = 0 and re- 
placing j by k +}. 


I am indebted to R. G‘. Stanton for pointing out the general form 
of the inversion formula used in Theorem 2. 
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POPULATION GROWTH UNDER THE INFLUENCE 
OF RANDOM DISPERSAL 


H. D. LANDAHL 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


A population is considered which grows according to the logistic equa- 
tion while spreading out atrandom. An approximate method is used to ob- 
tain transient and steady-state values for various simple boundary condi- 
tions such as that of a population started in an infinite one- or two- 
dimensional region with or without reflecting or absorbing barriers. An 
approximate steady-state solution is given for the one-dimensional case 
of two neighboring regions having different growth rates, mobilities, and 
degrees of attractiveness, 


In the present paper approximate expressions relating the varia- 
tion of the density of a population with time and place are derived 
for the case in which the growth rate is of second degree in the 
population density and the population spreads out by a random 
process. Such a situation has been treated by Skellam (1951) for 
certain cases in the steady state. 

Let n(p,¢) be the number density of the population under consid- 
ation, p being the distance from an origin and ¢ being the time. In 
two-dimensional cases we shall restrict ourselves to central sym- 
metry. Let a and f be the coefficients of the linear and quadratic 
terms in the expression for the growth rate as a function of n. De- 
note the diffusion coefficient by D = y?/2vT, where y is a mean 
‘‘jump’’ distance, 7 the time between jumps, and vy the number of 
dimensions of the region involved. Then if D is a constant, n is 
given as a function of position and time by the solution of the fol- 
lowing equation, together with boundary conditions: 


dn 2 2 

=——r(in—Bn +DV Ns 

Sri ee (1) 
1 y? ld 

Te an ot Ee 
QvT’ ar 
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n* being the value of n when the birth rate minus the death rate is 
zero. 

We shall first suppose that the initial distribution can be satis- 
factorily approximated by an expression which is linear in p and 
which has its largest value at p = 0. At some distance p =7), the 
population density is negligible so that we may put n(r,,0) =0. 
Denote n(0,0) by n). 

It is most likely that we would be interested in a population 
which is introduced into a fairly small region, the initial population 
density being relatively low. Hence we would generally be inter- 
ested in cases in which Bn? << an,. However, we shall limit our- 
selves only to the situation in which Bn? < an). 


Two-dimensional case for t<t*. Consider first the case in 
which the population can spread out over a two-dimensional sur- 
face. As long as an < Bn”, we may, to a first approximation, fol- 
low n(p,¢) by making n a linear function of p in which n becomes 
zero at some distance p = r(f¢). 


p 
n(p, t) = (0,1) fl ak (2) 
The quantities n(0,¢) and r(¢) remain to be determined from con- 
siderations of continuity (Landahl, 1953) in the following manner. 
Let n(p,¢) and n(p,¢+ At) be two successive distributions. Let 
M, be the total number of individuals at time ¢ which are within 
p = 1/2, while M‘, is the corresponding number at time ¢+ A¢. Let 
M, be equal to the total number of individuals which are between 
r/2 and r at time ¢ while M‘ is the number between 7/2 and r(¢ + A 2) 
at the time ¢+ At. Now M‘+M,-M,-M, represents the total 
addition to the population in the time interval A¢. This must be 
equal to the space integral of an(p,t)A¢—Bn?(p,t)At. In this 
way we find the following relation: 


d(nr*) 
dt 


3 
= nr? See o 37 = n(0,?2). (3) 


Similarly, if AM, is the increase in M, due to the space integral 
of (an — Bn®) At, while AM, is the corresponding increase in M,,; 
then M’.-M, - AM, or AM, — M‘,+M, is the increase in M, which 
results from the diffusion flow given by Dn(0, 2) divided by r and 
multiplied by the area. A linear combination of the resulting rela- 
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tion and of expression (3) leads to 


5 dr 8D Le a A 
het ae B nr. (A) 
If we make D =0 in (4), then 740, whereas actually n must 
eventually become equal to «4/8 everywhere within the region fi. 
but remain zero outside. The linear approximation is evidently un- 
satisfactory. A better approximation is given by a function which 
is constant for a certain distance, then decreases to zero at re 
We shall limit ourselves to cases in which the last term on the 
right-hand side of (4) is less than the first, so that to a first ap- 
proximation we may drop the last term, integrate, and obtain 


r= Vr? +8 Dt. (5) 


If N(t) is the space integral of n(p,¢), i.e., wr?” or M, +My, 
then (3) together with (5) leads to 


at 
Nae 


N(é) = 6 
(2) 3N,e74 7 = “yp (6) 
apr ets + at)- 1 
where A and X are defined by 
re ao, Tp is Th? 
LSS AeTe N, = 71,°7(0) = ; n(0, 0); (7) 


and Ki z is an exponential integral defined by 


x 
Ei o- f we dv. 
ie) 


Expression (6) can hold only if n(0,t) <%/B =n", the latter 
value being that value at which an =£n?. Let ¢* be the time at 
which the equality is reached, i.e., when n(0,¢*) =n". Then we 
find from (6) that for ¢ < ¢* 


3N,e% (1 - p/Vr,” + 8 Dt) ‘ (8) 


-A 
SB Noe" pi (4 4 a) - Bi Al} 
16 7D 


n(p,t) = 
7(t,7 + 8D?) fi + 
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Making n(0, ¢*) = n*, we have an implicit expression for ¢*. We 
shall introduce the following notation: 


* 


n 
7 alte a pe SKE) Se 
Then equation (8) may be written thus: 
ga ek 1 eA 
—--—KEi —Ki A- FR —=0. 10 
oye ee 7 (10) 


The graph of « ¢* obtained from (10) as a function of A with R as a 
parameter is given in Figure 1. If R > 1, a 2¢* can always be ob- 
tained from the figure. But if P <1, i.e., if the population density 
at the origin is initially greater than the natural limit, then « ¢* 
can be obtained only if A is less than a certain value which de- 
pends upon F (dotted lines, Fig. 1). Equation (10) has no real 
positive roots for «¢* when R<1 and A>2. For this case, the 
approximation which leads to (6) fails from the start. We consider 
only values of R and A for which « ¢* is real and positive. 

When & <1 and when A is less than a certain value which de- 
pends on Ff, there are two real positive solutions. The smaller 
represents the time when n(0,¢) crosses n* from above (broken 


FIGURE 1. 
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parts of curves in Fig. 1). The second and larger solution is the 
one we need. 


It may be noted that over a very wide range of values, o ¢* can 
be approximated by 


az“ = 2.5 (1+ log,, R/A), (11) 
provided that A < 5and R2A. Similarly we have 
r= r(t*) + 0.9 LV/1+0.4A+ log, RA, (12) 


where L, a characteristic length to be used later, is given by 


6D = 
pus V == 2¥3 = (12a) 


One-dimensional case for t<t*. If we are interested in the 
spread of a population in a long, narrow region, the results obtained 
from a one-dimensional approximation may be useful. This case may 
correspond to the situation in which there is a reflecting barrier at 
the origin. Or there may be no barrier but instead there is sym- 
metry with respect to the origin. For convenience we shall con- 
sider only half of the region. Following the above method we find, 
instead of (3) and (4), the following relations: 


d(7r) oS GA £, 
Ms = amr ~~ Bn ry 2n=n(0, 2), (13) 
2Qr dr 2 
= + 3 B (14) 


Using the same restrictions as those involved in obtaining (5) we 


find, for ¢ < ¢*, 
N, (1 - p/V'19 + 8Dt) e* 


SEIS Fee aye ae ee 
Vr? +8Dt jp + peter e*a~2dau 
Q 3 V8Da J, 


n(p,t) = (15) 


where N, = 7(0,0)7,, the initial number. The total number N (2) at 
any time ¢ < ¢* can be obtained from 


N(t) =3n(0, 2) Vr? + 8Dt, (16) 


and t* can be obtained from 7, where 7 is given by 


” x ~ 
e 2 pne dz ke Oh (17) 


Ga 
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The solution of this equation is facilitated by the substitution 
y=Vvn- 


It can be shown that for y 2 2 we have 


2 
ne*~ dz 0 aaa 1 e” 
wf dem 1) Ao 
0 Va 0 y 
where 
1 1-3 ino o 1.3.5.7 
Hdy) 2) te eee e 


—____ + ———____ 4 ——_______ + 
Dy?” (2y2)? " (2y2)8Q! (2ya) 43! 


Some values of f(y) are given in the following table: 


AG) 1.023 1.010 1.003 1.000 
f)H(y) .610 557 537 523 -5172 .5107 -5025 


For values of y < 2, the integral is tabulated by E. Jahnke and 
F. Emde (1945, p. 32). If we expand the exponential in the above 
integrals and integrate term by term, we find that the resulting 
series converges rapidly for y < 1. 

Using for y < 2 the values from Jahnke and Emde, and for y > 2 
the values which are calculated from the above table, we find from 
equation (17) a numerical relation between y and A, with R as a 
parameter. Because of (9) we thus also obtain the values of ¢*. 
The results are shown graphically in Figure 2. It can be seen that 
over a very wide range of values of # and A, including values such 
that A < 3, R > 1, we may approximate « ¢* by the expression 

R?2 
az" = 2.14 1.2 log,, = (18) 


One- and two-dimensional cases fort>t*. After the time ¢* we 
may attempt to approximate the solution by a function which is con- 
stant up to a certain distance s and decreases linearly between s and 
r. The function is illustrated in Figure 3 and defined for ¢ > ¢* by 


n(p,th=n*, O<p<s; 
n* (7 ~ p) 
(r=)? 
n(p,t)=0, p2r 


n(p,t) = s<p<nr (19) 
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R=j 
R=] 
R=] 
Yo Ve 
001 01 1 1 10 100 
A- 
FIGURE 2, 


The quantity s, like 7, is a function of ¢. At the time ¢+ AZ, the 
function will be of the same form but r and gs will have changed to 
7+ Ar and s+As. The area between the two curves represents 
the gain in the population, which is also the gain in the space in- 
tegral of (an - Bn?) At dp. We shall see subsequently that the 
two-dimensional case can be represented quite satisfactorily by 
the solution of the one-dimensional case for ¢ > ¢t*. We therefore 
consider the one-dimensional case first. In this case, equating 
the above two increases in population during AZ leads to an ex- 


s(t) S+AS 


a(pit)> —> 


5} s(t) r(t) r(t+at) 
p— 


FIGURE 3. 
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pression which, on dividing through by n™ and by Aé, may be writ- 
ten as follows: 


== 8) =F (r- 9), (20) 


where the dots denote derivatives with respect to time. 

We imagine, for a moment, that multiplication stops for a short 
interval of time A¢. During this interval of time n(p, 7) is deter- 
mined by diffusion processes alone. In that case n(p, ¢) is approx- 
imated by the dotted line 0,0, in Figure 3. The areas 0,0,0, 
and 0,0,0, are in this case equal because of the constancy of 
the total number. On the other hand we may consider that during 
the interval Az the diffusion stops while the multiplication con- 
tinues. In this latter case the line 0,0, will rotate clockwise 
about O, to 0,0, because of the increase in the total population. 
It is seen that the end result of the two situations 0,0, is not 
changed, if both cases occur simultaneously, provided that At — 0. 

The area 0,0,0, represents the total number which moved to 
the right of O, in the time A¢ due to diffusion. But this number is 
also given by the product of the diffusion coefficient times the 
gradient of the number density times the time A¢. Equating these 
quantities and dividing through by n* and Aé, we find 


ae D 


(21) 
Introducing 
Y= r—s, (22) 
multiplying (21) through by 4 and subtracting the corresponding 
terms from (20) we find, on rearranging, 


d(x?) 2 
tt == a(L? - 2%), (23) 


where L is defined by expression (12a). 
Since r = r* for ¢ = ¢*, the solution of (23) is given by 


. - 22) 
=(r—s)?=L2-(r*2-L%)e 3 : (24) 
Substituting (24) into (21) and integrating, we find 


tee tog fe ee ee , 
on ie te + L) ( [7,2 “ (r*2 = L*) e72 Ut t*)73 _ 


“* 


(25) 
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With r(¢) thus given, s(t) can be easily obtained from (24). With 
r(t) and s(¢) calculated, the approximate value of the population 
density n is given by (19) for any time ¢ and distance p- Note that 
after a sufficiently long time 2=r- 8 approaches L and f ap- 
proaches ¥20D/3, so that the population wave has a constant 
shape and velocity. 

We next determine under what conditions the above results can 
be used for the two-dimensional case. Using the same argument 
which leads to (20), we find for the two-dimensional case 


Qer+re+ sr+2ss=1a(r—s) (rt 8). (26) 


Proceeding as above we now find instead of (21) an expression 
similar to (21) except that it has the following factor on the left- 
hand side: 

(r+ 8)? +1rs+rs 


(r+ 8)2 + 73+ 3? Eye 


This factor equals unity for r= s and for s = 0. Now s is less than 
r. When s = 0.297, the above factor has its maximum value, 1.10. 
At this time, 7 increases at a rate which is ten per cent slower 
than that given for the one-dimensional case. Therefore we shall 
neglect this error and introduce (21) and (22) into (26), eliminating 
s and r. We then find an expression which is similar to (23) except 
that the right-hand side is multiplied by a factor 
1 

3 : (26b) 
8 
The limiting value of a for ¢ = ~ is again L, but at ¢ = ¢* the deriv- 
ative of 2? increases or decreases at a rate 14 times as fast as 
that of the one-dimensional situation though this excess rate soon 
disappears. Since this factor affects only the rate of approach of 
a to L, and this is most important near the beginning, we may in- 
troduce a numerical coefficient between 1 and 1.5 on the right-hand 
side of (23). This in effect amounts to multiplying a by the se- 
lected factor in the subsequent expressions. When #(¢*), the ini- 
tial value of z, which is equal to r*, is fairly close to L then itis 
clear from (23) that # remains almost unchanged and the factor be- 
tween 1 and 12 is unimportant. 

From expression (12) it can be seen that r* has a value greater 
than or about equal to 0.9 L, this value increasing very slowly 


T— 
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with R/A. The ratio A/F is proportional to the initial total popu- 
lation number N(0) and is equal to 9N(0)/N,,-where N, is the 
total number of individuals that are contained in a circle of radius 
L when n is everywhere equal to n*. 


One-dimensional case, central symmetry, absorbing barrier. Let 
there be an absorbing barrier at a point Z and a reflecting barrier 
at the origin, or let there also be an absorbing barrier at —Z, the 
initial distribution being approximated by a triangular one, sym- 
metric about the origin. Consider first the steady state. In this 
case n(Z,) is zero since L is an absorbing barrier. Suppose that 
n(p,e) is a linear decreasing function of p, which is zero at p = Z. 
Then in the steady state the integral of an — Bn? over p, which is 
the net increase per unit time, must equal the outflow Dn(0,~)/Z 
due to diffusion. Equating these expressions we find 


; AE an 
n(0,~) = n* (~- =) . (27) 


For Z < X = /2D/ua, n(0,~) is negative, i.e., there is no solution 
except n(0,0) =0. For Z>/3X =\/6D/q, we have n(0,~) > 1. 
Hence there must have been some time ¢*, at which r(¢*) = r* < \/3X, 
when the linear approximation should have been discarded. If this 
is done, we find that the steady-state solutions are as shown in 
Figure 4. When Z < X, the population becomes extinct. For Z > X, 
the population density can be approximated by the functions shown 
in the figure, the position of the barrier (Z) determining the shape. 


1.0 


1X V3x 2x 3X 
2 —_ 


FIGURE 4, 
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If Z > /3X, then n(p, 0) = n* for x < Z—\/3X. Note that the only 
parameter is (%/D). In this situation we may compare the value of 
X with that of the exact solution for the steady state. It is equal 
to (7/2\/2)X = 1.11 X in our notation (cf. Skellam, loc. cit., p. 
204). In Figure 4, two exact solutions are included (broken curves) 
to illustrate the nature of the approximation. 

The transient solution will depend on the value of Z. We sup- 
pose that 7,, the range of the initial population, is less than Z. 
Then 7, the range at the time ¢, will increase according to (5) until 
either r(¢,) = Z or n(0, ¢*) = 1, depending on whether ¢, is smaller 
or larger than ¢*. During this time n(p, 2) is given Ey expression 
(15). Ifz, <2", fen after ¢, when r= Z, we find the expression 
for the average number density 7 = +n(0,t) in the following way. 
Subtracting the diffusion flow Dn(0,t)/Z from the integral of 
an —n? and setting this equal to the rate of increase of 7, we 
obtain 


i -— Bn. (28) 


Note that the steady-state expression (27) follows from expression 
(28). If Z< xX, then it is evident from (27) that n decreases to 
zero. The decrease is a simple exponential if the last term is 
negligible. Using the boundary condition 7(t,) = $7(0,¢,), where 
t, is given by Z=\/r? + 8D¢, and n(0,¢,) is given by (15), we 
find that the integral of (28) is the desired solution. IfX¥ <Z< V/3X, 
7 increases but the above indicated solution still holds. If Z > /3X, 
then the above~indicated solution is satisfactory until a time ¢, 
when 27(¢,) =n". Thereafter a solution is obtained by assuming 
that n is equal to n* in the range from p = 0 to p = 8, and that n 
decreases linearly to zero from p= 8 top =Z. The quantity s is 
then determined from the requirement that the space integral of the 
growth rate on — Bn?, less the diffusion flow Dn*/(Z — s), equals 
the rate of increase of the total population. 

If the time ¢, is greater than t* obtained from (17), then from 0 
to ¢* the solution is given by (15) and from ¢* to ¢, by (24), ¢, be- 
ing the time at which r= Z. Thereafter the Bolton can be Fond 
as indicated at the end of the above paragraph. 


One-dimensional case, central symmetry, reflecting barrier at Z. 
We consider next a case like the previous one but in which there is 
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a reflecting boundary, e.g., a physical boundary which cannot be 
passed. If the distance Z to the boundary is less than r*, the 
solution is the same as if the boundary were absent as long as 
t<t,<¢*. After ¢, let the distribution be given by a linear ap- 
proximation in which n = €4/B atp=0 andn=nn* atp=Z. The 
quantities € and 7 can be determined as functions of time from the 
following conditions. The increase in population in the region from 
0 to Z/2 is represented by the integral of the net birth rate from 0 
to Z/2 less the diffusion flow out of the region. The net increase 
in the region from Z/2 to Z is represented by the integral of the 
net birth rate from Z/2 to Z plus the diffusion influx from the first 
region. Putting 


w= E+, V= — hs (29) 
we find 
. 2S L? 
des © heat eae eg I (30) 
ae gy ree 1 1 
(= ——— aS arn a 7 
5 ( 5 = 0) (31) 


At t=t,, we have n = 0 and = € = n(0,¢,), where n is given by 
(15) and ¢, = 3(Z? -1,°)/aL?. Hence att=¢,,u=v=&. With 
these boundary conditions the solution of (30) and (31) for u and 2, 
and hence for € and n, gives the desired approximate solution. If, 
in addition to the requirement Z < r*, where r* is of the order of 
magnitude of L, we also require that L? be substantially greater 
than 3 Z? or, say, Z < L/3; then from (4) we find 


v= €, eh AZ, (32) 


where ¢’=¢~t,. Now the last term in the parentheses on the 
right side of (31) is at most €,/6, relative to u, and &, is less than 
1, It is clear from (31) that u approaches 2, i.e., = =1. The 
last term can have an important effect only when uw approaches 2, 
in which case u~-iu? becomes small. But then n will already be 
fairly close to its steady-state value. Therefore, we shall neglect 


this term, Solving for w from (31) and using (32) and (29) we find 
for Z<4L,Z <1", 


yf bes at” é 
ae PM eee ee eee 33 
arr f + (2- &) err aa 
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n= = 
by + (2£,) 672 o7 
IZ > r*, so that ¢, > ¢*, then the solution for ¢ < ¢* is given by 
expression (15), while for ¢* <¢< ¢, the solution can be obtained 
from (25). After ¢, a very rough approximation can be easily ob- 
tained by assuming a linear distribution which is equal to n* at 
s*, corresponding to r*, and to some value 7n* at p=Z. The 


value of 7 depends on the time ¢. Since s* is not permitted to 
change, we find the following expression for determining 7 


n= (1-7) (1+ 27). (35) 


Since s* is given by (24), (25), and (22), and n = 0 at ¢ = ¢, there- 
fore if ¢’=¢-¢,, 7 is given by 
ett 4 


| BB remer mana one 3 
ng 


1 iis, 2 
mae Oat es (34) 


(38) 


With this value of n the approximate solution for n is then known 
for ¢ > ¢,. 

The above solution is a linear approximation and therefore shows 
a negative slope at the reflecting barrier. Actually the derivative 
must be zero. Similarly at the origin the derivative of the number 
density must be zero. The present method attempts to give only 
the best linear approximation to the number density at any time. 


One-dimensional case with two regions having different param- 
eters. Consider a linear region in which to the left of a point Y 
the parameters are D’, «’, 8’, while to the right the parameters 
are D’’, «’’, B’’. Let a population be introduced far to the left so 
that after a time there is a wave of population growth moving to 
the right with a velocity v’ = /2’D’/3, the growing region being 
of width L’. When the front reaches the origin new approximations 
must be introduced. But after a time there will be a wave of width 
L’’ moving along the axis to the right with a velocity o’. Con- 
sider the population distribution in the neighborhood of the origin 
when the growing region is far enough to the right. Define F by 
the relation 

peal 


: ait (37) 
a 
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Then we may approximate the steady-state solution as shown in 
Figure 5 (heavy line), where F is chosen to be positive. The 
quantities s’, s’’, and f are determined by the following three con- 
ditions. (1) The flow from the right across a line at Y must equal 
the flow to the left, i.e., D ,f/s’ = D,(F - f)/s°’’. (2) The integral 
of the excess growth between Y and s’” must be equal to the out- 
flow to_the left. (3) The deficit of the net growth in the region 
from — 8° to Y must be equal to the flow in from the right. Without 
appreciable complications we may introduce a partition coefficient 
y between the regions. This coefficient measures the preference 
for one region over the other. Such a case may occur where a 
wooded area borders a meadow. In this case the quantity n/n* will 
be discontinuous at Y, being (1+ /) on the left and y(1+/) on 
the right. These conditions lead to the following relations: 


Air #1 
ee (38) 
yAo +a 
; : (39) 
oO ————T, 
v1+if 


ia V 3(1+ F) 
oO = — 
14 2yeF 4 ops, “ 
c 3” es (hy Ct 
oO ~ ye? Co “yer? A-= sates (41) 


, 


where 


The values of f, 0’, and o”’ can be obtained by successive approx- 
imation as follows. Let o)’ be the zeroth approximation to o’ 
which is obtained from (39) by setting f = + F. Similarly let o)’’= 
o°’ for f= 4F. Then /,, the first approximation for f, is given by 


FIGURE 5, 
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(38) with o” = ia and oe =o). Then o,’ ando,’’, the first ap- 
proximations to o’ and o*”, are given by (39) and (40) with f = ie 
Successive approximations converge rapidly unless the ero regions 
have very different parameters. The values of s’ and s’”’ are not 
greatly different from X’ and X’’ respectively unless the regions 
are very different in regard to the population density that can be 
supported. These values are therefore less than L’ and L’’ (L = 
V12X). Hence, the front of the population wave reaches the ori- 
gin (¢,, Fig. 5) before the crest (where n becomes equal to n* for 
the first time) reaches the point Y - s’, After the time ¢, we may 
approximate the solution as indicated in Figure 5, a separate solu- 
tion being required for each interval between ¢, and ¢,- In each 
case it is necessary to divide the regions into about equal parts 
to determine the flow from one half to the other. 

If in the above example we find that o’ X’ > Y, then it is nec- 
essary to suppose that even at the origin, which is a center of 
symmetry or a reflecting boundary, the population density exceeds 
n’* by some amount 6n’*. Following a procedure similar to that 
discussed above, we now find an expression which is similar to 
(38), except that o’ is replaced by Y/X’ and that 60°’ is added to 
the numerator on the right-hand side. Instead of (38) we find the 
following expression involving 6: 


Y7(f(1+2/f)+20+ 36? +f 9] = (f — 6) X*?. 


Since generally f will be considerably larger than 0, the quadratic 
term in @ will not appreciably increase the number of steps in the 
successive approximations for the estimation of f, 6 and o’’. The 
expression for o’” is the same as in (40). 

If we consider the two-dimensional case with central symmetry 
where the boundary between the two regions is at a distance Y 
from the center, then if Y > Y*, Y* ~ X’, we find the approximate 
solution from three equations, the first of which is identical with 
(38). The second equation is similar to (39) but has 1 +2f- 
o X’(1+4/)/3Y instead of 1 + 2 f under the radical sign. eee 
pg <4 ar Y & X’, the correction term is not likely to be impor- 
tant, unless f is large. The third expression is similar to (40) 
but with o” X’(1+F +y+yf)/6Y added in the denominator of 
the radical. As long as Y is not too much less than 8°’, this addi- 
tional term does not appreciably complicate the calculation of 
o’,o’’, and f by successive approximations. 
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If data are available and if they can be approximated by four 
linear parts as in Figure 4 (though there may be a discontinuity at 
Y if y #1) then, since F, y, f, 3’, and 3°’. are given directly, we 
may use expressions (38)-(41) to calculate X’, X’’, and A. From 
these we directly obtain D’*/D’. and «**/a’. for either the one= or 
two-dimensional case. If any D or & is independently known, then 
the others are known. Thus a single additional point on the tran- 
sient process which involves either of the regions is adequate to 
determine all the parameters. 


This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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The logical network of a Pitts-McCulloch type is developed further 
with the specific addition of a motivational system which acts as a dif- 
ferential reinforcer to the network. A consideration is given to the vari- 
ous ways that probabilities enter the logical net transforming it into a 
probability net, and lastly some notes are added on the possibility of in- 
cluding the strategies of Theory of Games in the logical network. 


1. Introduction. A basic system of counters for an inductive 
logical network has already been described (George, 1956). This 
system was shown to be perfectly general and to be sufficient to 
count the conjunction or disjunction of any number of events what- 
ever, with respect to any number of inputs whatever. Any number 
of counters could be used. 

We now wish to add the mathematical equations for a motiva- 
tional system. At the same time we shall consider the possibilities 
of an alternate notation for these networks, and also consider the 
various ways in which it is necessary to introduce probabilistic 
notions into the system. 

The set of three conjunction and three disjunction counters, as 
part of a B(3,3) network (George, 1956), is reproduced here in 
slightly improved form. Thus Figure 1 shows the network where 
the letter B was chosen to represent a belief which had a special 
significance in a previously developed system of behavior (George 
and Handlon, 1955). 

The following are the conventions on which the networks are 
drawn. A circle represents an element which is in one of two states 
at any instant; it is either ive or dead. If live, then its output 
fiber will be live and sending out an impulse. The output fiber 
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must be in one state at any instant which is assumed to be the 
firing time of all the elements. Thus the output may bifurcate as 
often as is needed, but it can only be in one state for all its bi- 
furcations. There can be as many inputs as we like, and these in- 
puts can be divisible into two classes—excitatory and inhibitory. 
The excitatory ending is represented by a filled-in triangle and 
the inhibitory by an open circle. A heavy dot is used to indicate 
that wires which cross over also have contact, except when three 
wires are involved. In this case, there will be no ambiguity since 
no loose endings are involved except for the input and output ele- 
ments. The number A (1 or 2) in the circle represents the threshold 
number. If e represents the number of live excitatory fibers at any 
instant and z the number of live inhibitors at any instant, then the 
condition for the element to fire is e27z2+h. 


Il. An alternative notation. An alternative network, for the net- 
work variously called the B-net or C-system, as shown in Figure 1 
is now described. First we shall state the equations of the net- 
work in Principia Mathematica notation, suffixed for time. Then 
we shall repeat the same equations in a type of notation usually 
associated with the name of Lucasiewicz. In any event the nota- 
tion is derived from the propositional calculus and the lower func- 
tional calculus, suffixed for time. 


ifs Hg Opa i (1) 
Es So Oa 643 (2) 
k} = 4.055 oa, = (3) 


op(kap er dt_,) w(ohig (Ed et-~ (bw apy (4) 
In deriving (4) use is made of the obvious simplifying condition 
Terie (ij v Se (5) 
Then the generalized equation for any number of c-counters is 
of = (ky erty) v opt! v ((Et)e™(%)» ~ (12_, v 12). ©) 


The generalized equation here makes use of (5) and of the ob- 
vious condition that results from the fact that, if c™*! fires, then 
c” must necessarily be firing. Similar equations can be derived for 
the d-counters, and a final condition on the ‘‘key’? element (ab) in 
Figure 1. The primed notation always indicates an output element, 
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FIGURE 1 


and the combined primed elements are the key elements. These 
key elements will be composed of every combination of all the 
primitive inputs a, b,...,n that the system contains. The final 
condition on the key element (ad)’ is 


(a), = (4,_,° 5,1) v(a,_,° Coy) Vv (5,_4° Geng) ? (7) 


where 


ie Se ee (8) 


The alternative notation of the Lucasiewicz type would write the 
same equations in the following form: 


Q11KNaib1; (9) 
Q11Ka1N01; (10) 
Qk1Kalb1; (11) 
Qci1LKk11Nd11K c211K & tc11NL 111121; (12) 
Ck1L1112; (13) 
QcniLLKkiie(n - 1) c(n + 1)1K & tenNL 111121; (14) 
Q (aby LLKa1b1kaiciik 61c11,; (15) 


CeliNdl1. (16) 
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The conventions of the notation are mostly obvious, but the fol- 
lowing is a list of transforming equations from the one notation to 


the other. 
Qab = df-a= b; 


Cab =df-a> b; 
Kab = df-a- 6b; (177) 
Na=df-~ a; 
Lab = df-av b; 


S tf(t) = df-(Et) f(t). 


Now we should notice also that cll means Fite c21 means Gout: 
and so on. Should any ambiguity arise from such notation, paren- 
theses could be used as exemplified in (14), where they are not 
absolutely necessary. Equations (1)-(8) are thus repeated in (9)- 
(16). It need hardly be added that, after appropriate substitution, 
we can state the condition of firing of the key element in terms of 
the input elements alone. This matter will not be pursued here; 
neither will that of the very simple generalization of equations for 
any number of inputs. Both these matters have been discussed 
elsewhere. 


Ill. Probabilistic considerations uithin the deterministic net- 
work. The design of the counting system already outlined is clearly 
the design of an inductive logical machine. Thus it follows that 
the basis of the network is probabilistic. Indeed the method of 
connecting the counters, for scoring and erasing purposes, will 
directly influence the sort of probability that will be realized by 
the network. A further influence on probability considerations will 
be the number of counters available, relative to the number of 
events to be counted. Let us consider the second point first, for 
a network of the type already outlined in this paper. 

It will be easy to see that such a machine, or network, as the 
one under discussion, even with many more counters, will remember 
only the recent past if the number of events it has to count is large 
relative to the number of counters. If the number of such events is 
relatively small, then a frequency probability of the Laplacian kind 
will be the outcome. We can state these results in two theorems. 


Theorem 1. If % 28 the number of events that the network is to 
count and n is the number of available counters, then if « <n/2, 
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the count will be on the basis of a Laplacian frequency probability. 

The proof of this theorem follows immediately from the nature of 
the connection of the counters, as shown by the equations that 
govern their firing conditions. Half of the n counters are assumed 
to be for counting conjunctions and half for disjunctions. Thus, 
if a <n/2, even if the whole count was either of conjunctions or 
disjunctions it would not exhaust the full set of counters and thus 
no events would occur that would not contribute to the final state 
of the machine which will fire its key element according to whether 
or not 


This clearly is the same condition as 


n m 
1 
Nas 5 he eh a 


i=1 k=1 


where > — i c, + Da d.; 
k =1 i=1 j=1 
and this is of the form p = m/n, where p is the ordinary Laplacian 
probability. 

It will be noticed that at the final state there will be no way of 
telling in what order the events occurred that contributed to that 
final state. 

Corollary 1. If « >n/2, then the conditions of theorem 1 will 
still be fulfilled if the number of conjunctions (or disjunctions) 
that make up the events is less than n/2. 

The proof of this is obvious, and the corollary merely imposes 
the condition that the number of conjunctions and disjunctions 
should be roughly the same. 

Corollary 2. If « >7n then, regardless of the ratio of conjunc- 
tions to disjunctions, the probability count will be weighted for 
recency, and some events will not influence the final state of the 
machine. 

We only have to consider the case of n conjunctions to realize 
that the n/2 conjunctions that came last will not have further 
changed the counting state. 
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Theorem 2. The maximum memory of the network is given by the 
number of conjunction-counters or disjunction counters, but not 
both, that are firing at any instant. 

Thus, if at some ¢ the network has r live counters (they must be 
either all c- or all d-counters), they might have been compiled in 
any way whatever. The only knowledge we have, without knowing 
the previous state at a previous time, is that there has been an ex- 
cess of one type of event over the other, assuming that the network 
started with zero count. In the last r instants there could have 
been n/2-—r disjunctions, if r> n/4, or there could have been r 
conjunctions, or any combination whatsoever of these, according 
to the state r instants ago. 

It is clear that by the memory we do not mean the most recent 
events that have affected the machine, but the excess of one kind 
of event over another. 

All that has been shown above makes it clear that the selection 
of the number of counters relative to the number of events to be 
counted can weigh the count (of frequency) in favour of recency 
as much as we want. 

We shall now consider the method of connecting the counters. 
What has been proved holds for the system as depicted in Figure 1. 
Naturally we would change the conditions by changing the counter 
connections. By suitably chosen counter connections we could 
realize any desired form of the probability function. We shall not 
pursue this matter in any detail, but we shall mention one simple 
modification in Figure 1 that would immediately change the nature 
of the counting. If we remove the connections from ¢! to d', and 
d' to c', which are inhibitors, then we should have the effect of 
stopping the practice of not counting disjunctions while conjunc- 
tion counts exist, and vice-versa. In such a new system, a dis- 
junction would be added and a conjunction deducted at the same 
time, or vice-versa. 

It is clear that by suitable choice of connections, coupled with 
suitable choice of the number of counters, we have a very wide 
scope in the possibilities of counting. In this connection, careful 
attention will have to be paid to various approaches to probability, 
induction, degree of confirmation, and degree of factual support 
(Hempel and Oppenheim, 1948; Kemeny and Oppenheim, 1952). 


IV. Multiplewing. The method of multiplexing that was suggested 
by von Neumann (1951) could be applied to the nets that we have 
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described. This would increase the probability of their working 
effectively, even though the basic components were not always 
working effectively. Again we do not intend to give a detailed 
analysis of a multiplexed system, but it is clear that we can fol- 
low the principle of the executive and restoring organ if we dupli- 
cate the network wherever it is possible. A detailed network could 
always be supplied that used the same majority organs as were 
suggested by von Neumann. 

In Figure 1 it is clear that a duplication of both the inputs from 
a and 5 to &* with an increase of the threshold to 3 would permit 
the failure of any one of the inputs to fire at any instant and would 
still bring about a correct count. This is an example which can be 
applied without much difficulty to all the other elements. Only the 
l-elements would be recalcitrant under such changes, although 
even these would then permit some malfunctioning and still be 
capable of keeping a proper record. 

Another matter that will receive no detailed consideration in 
this paperis the use of random elements. These can be used almost 
anywhere at any time, and the predictability of the network will 
largely be destroyed as a result. Their use in explicit design may 
be somewhat limited, even when the network is explicitly intended 
as a model of organic behavior, since the fluctuating interference 
of the emotional systems on the cognitive system is not a random 
matter, but a systematic causal one. Random elements would at 
best be a makeshift device for which a deterministic substitute 
should ultimately be found. The indeterminacy of the system would 
be sufficiently supplied by the interference between one system 
and another, and the attempt to retain proper communication would 
be made in terms of multiplexing. 


V. Motivation. Any cognitive system of counters that we have de- 
scribed would systematically count everything that occurred which 
it was capable of discriminating, unless there was some system for 
selectively reinforcing certain events at the expense of others. 
This certainly happens in organic behavior. Motivation, which is 
basic to all organisms, depends primarily on survival and immedi- 
ately on the preservation of organic homeostasis by eating, drink- 
ing, etc. 

A certain subset of stimulus-response activities is assumed to 
be innately satisfying, or unsatisfying, and is built-in. All the 
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other stimulus-response connections gain motivational content by 
association with those that are built-in. 

We may now give a general description of the motivational sys- 
tem. Figure 2 shows a particular M-system with just two m- 
chains composed of 4 elements per chain: m', m”, m®, m*, n', n?, 
n® and n*. 

In-general we shall say that a motivational system is composed 
of n chains of m elements each. Let these be written: 


m 


mm 
x m*, = ge etc., 


i=1 j=1 


which may be written M,N,... for short. 

Then an m-chain M, say, fires if and only if more than some 
number of elements become live. Let this number be a fixed in- 
teger k. (It is clear that we may want to make this a more flexible 
condition dependent, say, on rate of ellicitation.) 

If r, the number of live motivational elements, is greater than k, 
then M fires into the cognitive, or Belief-system. The nature of the 
connections that the M-system makes with the cognitive system 
will be discussed briefly a little later. 

Figure 2 shows that delay elements p', p?, p®,... will be neces- 
sary, suitably chosen according to the delay expected between the 


To FROM 
C-SYSTEM C-SYSTEM 


P3 


FIGURE 2 
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response and the recording of the result of the response. The de- 
lay elements allow the result in terms of the M-system to be con- 
nected with the counters for the appropriate responses. This 
process of selection makes necessary the selector elements drawn 
in two columns. These elements s', s?, s%,... are arranged in 
pairs so that two such elements are associated with each possible 
response. It should also be mentioned that there can be many more 
delay elements than are strictly necessary to associate a particu- 
lar response with a particular outcome. This has the effect of re- 
inforcing other responses (or stimuli) that are closely associated 
in time with the ones that are being reinforced, because of the 
satisfactory outcome of the response. The actual choice for a 
machine would normally allow some leeway in this matter and make 
such a gradient of reinforcement inevitable. 

The selector elements can be divided into the left-hand and 
right-hand members where, as shown in Figure 2, the left-hand 
element is the one that will fire if and only if the m-chain is not 
firing; the right-hand element fires if and only if the m-chain zs 
firing. 

The general condition for an m-chain to fire is 


, = (Et)m(d), (18) 


where m(Z) = >» m, and r>k, k being the constant referred to 
r=1 

above. Otherwise M does not fire at ¢. Then if some a fires so 

that a’, (ab)’, (ac)’,... etc., fire, 


(8.4 806i 'Saci’ °° 318.4 Cr errr eee = (19) 


where A means the exclusive sense of ‘‘or,’’ and the suffices 1 and 
2 mean ‘‘left-hand’’ and ‘‘right-hand’’ respectively. The letters 
a, ab, ac,... refer to the associated responses. 

The principle of the M-system is simple. A condition of a 
logical type, as those in (1)-(16), could be given for every element 
in the system. This will not be attempted here, but it should be 
clear, from what has been said, just how these equations should 
be written for any particular case. Figure 3 shows an alternative 
M-system with loop elements for m- and n-elements. 

Now to return to the linkage between the M-system and the C- 
system. There are many equivalent ways in which this can be 
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FIGURE 3 


done. Perhaps the simplest is to dispense with the £- and /-elements 
and join the wires from the left-hand selectors directly to the. con- 
junction counters, and the wires from the right-hand selector ele- 
ments directly to the disjunction counters for excitation, and the 
other way around for inhibition. A condition on the firing of the 
delay element p' of (9) can be made by connecting the m-chains 
directly to this delay element. 

We shall not give the full set of equations for the new system, 
but we give the equation for c' alone, and all the others can be 
derived in similar fashion. Thus the generalized equation for c? is 


QclLLLKK(LL...£8118381...)N(LL ...1821841...) (20) 
Nd11K c21% tce1tkk (LL ...Ls11881...)c21% écll, 
where we have used s‘, s*,... for the left-hand elements and s?, 
s‘,... for the right-hand elements. 

The motivational system will not be developed any further here, 
but the question of the neurophysiological interpretation is of the 
the first importance and should be undertaken in the future. 


VI. Theory of games for finite automata. The last note we should 
add is a general one on the strategies suggested by Theory of 
Games for our logical networks or finite automata. 
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Consider the utilities matrix 


11 in 
Fo Qo Gon 
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of utilities a;, which measures the utility of row ¢ for column j. 
The rows represent possible strategies from which the individual 
chooses and the columns represent states of nature. The return to 
a player choosing row z when nature is in state j is a;;- Then 
there are various rules supplied for the strategy to be played 
(Thrall, Coombs, and Davis, 1954). 

Four of the best known such strategies are those of Wald, Sav- 
age, Hurwicz, and Laplace. Wald’s principle of minimax states 
the principle as maximizing the minimum a; ; Savage states the 
minimax regret principle by defining the regret matrix ras 
max @,;. The matrix r;; measures the difference between what is 
actually obtained from a choice and what could have been obtained 
if the state of nature had been known. Wald’s principle is then 
applied to the regret matrix. Hurwicz’s system depends on the 
choice of a coefficient of optimism 2, say; for each of the rows let 
A be the largest and a the smallest component, and then maximize 
the function zA + (1 - 2) a. 

The last and best known strategy is that of Laplace and in- 
volves the choice of the row with the highest average. The fol- 
lowing matrix due to Milnor (see Thrall, Coombs, and Davis, loc. 
cit.) illustrates the differences. 


2 oben OP cel Laplace’s choice 
bortinele «1 Wald’s choice 

Ogrtinith> 0 Hurwicz (x = +)’s choice 
Lb whiad Savage’s choice 


There are many things that should now be said to develop stra- 
tegies from these simple beginnings in game theory, but here our 
aim is merely to indicate the link between game theory and logical 
networks. 

First we should notice that these strategies are related to sub- 
jective probability estimates and guessing behavior and apply to 
behavior at the level of uncertain information. This means that it 
will be of special importance at the initial stages of learning. 
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It is already clear that behavior as a whole tends to lend itself 
to treatment as a Markoff process, since our logical nets make be- 
havior a series of S-R-S activities. Those activities have proba- 
bility relations of the backward looking variety, which characterize 
a Markoff process, with respect to nature as an ergodic source. 

The behavioral question is which strategy will be employed 
under various conditions. The answer would appear to be that we 
may work out a best strategy on well-founded empirical evidence, 
but we want to see how the changes might be expected to occur in 
practice. 

The occurrence of beliefs of roughly equal value (George, 1956; 
George and Handlon, 1955) will give rise to states of uncertainty. 
Where, however, equality of beliefs exist—that is, an equality of 
strength in terms of some motivational principle of maximum reward- 
minimum effort—there will still remain the component variables 
that will generally differ. Thus a particular individual will evalu- 
ate such a state of affairs in terms of giving one sort of variable 
priority over another. Thus we are led to the pessimist’s and the 
optimist’s strategy and also other variations that reflect the ex- 
perience of the individual. 

In short, strategy techniques represent individual differences 
between different logical nets. They are the individual constants 
in behavior and are represented as indices or coefficients of one 
kind or another that represent such things as differences in per- 
sonality. They may be brought about in nets all of which start in 
the same manner, as a result of a different segment of environment 
being encountered or due to differences in the stability character- 
istics of the net, or both. By the stability characteristics we mean 
the extent to which ‘‘rapid counting’’ takes place when environ- 
mental conditions change a reinforced to a weakened connection. 
The motivational system hastens the slow counting change in the 
C-system and makes it possible to change habits quickly. 

In psychology the random activity of a Thorndike cat is followed 
by a definite strategy in escape as a measure of success enters 
the picture. If we change the rewarding conditions slightly, we 
may change the connections to some extent which makes the dif- 
ference between the use of one technique rather than the use of 
another one in its place. 

The problem, which now has to be described in detail for any 
particular net, is how that net acquires the strategies it performs, 
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relative to a particular environment. No attempt will be made here 
to carry through such a detailed analysis. 
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A NOTE ON THE GEOMETRIZATION OF BIOLOGY 


N. RASHEVSKY 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


In connection with a recently proposed topological model of an organism 
(Bull. Math, Biophysics, 18, 31-56, 1956), it is pointed out that the de- 
pendence of some properties of the studied topological spaces on the 
type of space in which they are imbedded reflects some aspects of the 
dependence of the organism on its environment, 


In a recent paper (Rashevsky, 1956) we suggested a possible 
*“‘geometrization’’ of biology, by studying geometrical, and in par- 
ticular topological, structures, the topological relations within 
which are logically isomorphic to the biological relations within 
an organism. We pointed out that one of the basic phenomena in 
biology is the process of selection of certain elements of the en- 
vironment, followed by subsequent assimilation. This selection 
by the organism occurs, however, only when the organism simul- 
taneously ‘‘loses’’ some of its elements either through catabolic 
breakdown which produces the energy necessary for the selective 
processes or by secretion of enzymes which act to a large extent 
as ‘‘selectors’’ on a molecular level. 

A selection logically implies the division of a set of elements 
in two classes: the selected ones and the rejected ones. It was 
pointed out that cases of division of a set in two classes by a 
geometric structure are quite common in mathematics. They are, 
however, insufficient for our purposes. What we need is a subset 
M of a set S, such that M divides S — M in two classes if and only 
if M loses some subsets of its own. Since, to the best of our knowl- 
edge, no such systems have been studied in topology, we invented, 
as an illustration only, a few cases of topological spaces which 
have the following characteristic: 
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If certain points of a subspace M of the space S are removed, 
then the points of S — M become divided in two classes. One class 
is characterized by the points acquiring some property P with re- 
spect to M; the other, by the absence of this property. A space M 
which is characterized by this peculiarity we shall call a A=-space. 

In the examples which we studied, the property P consisted 
either-in making M connected, or in reconstructing M by being 
added to the lost elements of M. 

This situation as has been discussed in detail in Joc. cit. pre- 
sents some basic logical analogies to the biological situations in 
organisms. It was also shown on an example how we may arrive 
at empirically verifiable biological statements by ‘‘translating”’ 
the geometrical properties into corresponding biological ones. 

The purpose of this note is to point out that the spaces which 
we have studied in loc. cit. have another characteristic which of- 
fers a formal logical analogy to organisms. Namely, the character- 
istic property of those spaces M depends on the structure of the 
space S of which they are subspaces. As an example consider the 
case discussed on pages 48-49 and illustrated on page 47 of loc. 
cit. The space M, considered as subspace of S = E?, shows the 
necessary characteristics. Upon removal of the points 8, and 6, 
which may lie on any of the subspaces A’, A’’, B’, or B’’ of M, 
any Jordan arc of S — M, which does not intersect M, which con- 
nects M’’ and M’”’ (see Joc. cit.), and whose end-points are cut 
points of either M’’ or M’’’, acquires the property of reconstructing 
a second space M, which is homeomorph to M but distinct from it. 

Consider, however, the same space M as a subspace of the 
space S, which is a strip in FE? contained between the lines pq 
and p’¢’ (Figure 1). Let M be situated with respect to S, as shown 
in Figure 1. Now, in order to achieve the same result, the points 
B, and 6, must necessarily lie either do¢A in one of the subspaces 
A’, A’’, B’, B’’, or they must lie both in either A’ B’ or both in 
A’ B’", If for example 8, lies in A’, while 8, lies in B’’, then 
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there is no Jordan arc in S,-M, which does not intersect M and 
which connects A’ and B’’, 

If we consider M as a subspace of S,, where S, is shown in Fig- 
ure 2, then M no longer possesses the characteristic studied above. 

Thus a given space # is or is not a A=-space, depending on the 
space S of which it is a subspace. 

All the above is perfectly obvious and to some extent trivial 
from a mathematical point of view. It is, however, much less so 
when we ‘‘translate’’ the geometric relations into biological ones. 

To the quality of being a A-space there corresponds the quality 
of an organism to select certain elements from its environment and 
to eventually reproduce itself. The space S—M corresponds to 
the environment of the organism. Translated into biological lan- 
guage, our above remarks amount to the statement that the ability 
of an organism to select elements of its environment depends on 
that environment. This at first also sounds quite trivial and ob- 
vious. The obviousness of this, however, is due to our familiarity 
with this fact which is so generally true. Of course, if the ele- 
ments which must be selected by the organism are absent from the 
environment, then zt is obvious that the organism cannot select 
them. But it is not obvious or self-evident at all that certain 
plants can select nitrogen from nitrogen compounds, but not di- 
rectly from air. If those plants are left in an environment which 
does not contain nitrogen salts, they will not grow in spite of the 
fact that nitrogen is present in air. 

If a general theory of A-spaces is developed, then the mathe- 
matical conditions which the space S (or S—M) must satisfy in 
order that a given M be a A-space, will give us, when translated 
into biological language, the biological or physico-chemical con- 
ditions which an environment must satisfy in order that a given 
organism may exist in it. But this will have definitely a predic- 
tive value. 


FIGURE 2. 
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REMARK ON AN INTERESTING PROBLEM 
IN TOPOLOGICAL BIOLOGY 


N. RASHEVSKY 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


The possibility of several homotopic classes of mappings of the graph 
of an organism onto the primordial graph (Bull. Math. Biophysics, 16, 
317-48, 1954) is considered. An application of this possibility is sug- 
gested for the theoretical determination as to what type of new biological 
functions may be acquired by certain cells which originally perform a 
different biological function. 


In our approach to the applications of topology to biology (Ra- 
shevsky, 1954; 1955a,b,c,d; 1956a,b,c) we attempted to represent 
an organism by a directed graph which is, so to say, the ‘‘organi- 
zation chart’’ of all its ‘‘biological functions.’’* The basic princi- 
ple of topological biology (Rashevsky, 1954) requires that the 
graphs which correspond to any higher organisms all be derived by 
means of the same multiparametric transformation from one or, at 
most, a very few ‘‘primordial’’ graphs which correspond to some 
primordial organism that either exists or existed in the past. Every 


*Although the semantic connotation of the words ‘‘biological function”’ 
as distinct from the word ‘‘function’’ used in mathematics was clearly 
stated in our first paper on the subject (Rashevsky, 1954, p. 330, foot- 
note), there have been misunderstandings due to the use of the word 
‘‘function’’ in two very different senses. A biological function is, for 
example, the process of digestion, vision, locomotion, multiplication, 
etc. The organism and its parts perform different biological functions of 
the above type. The biologist would say the function of the stomach is 
digestion, meaning that the stomach performs digestion. When we use in 
our publications the word ‘‘function,’? without the qualifying adjective 
‘*biological,’? we understand the word in the sense as it is understood in 


mathematics. 
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derived or transformed graph can be mapped in a many-to-one way 
onto the primordial graph. 

We have studied a number of different possible transformations 
and from them derived different biological consequences, some of 
which represent known empirical biological laws, others suggest 
new relations hitherto not studied experimentally. 

In all the above cases it has been shown how the transformed 
graph maps onto the primordial. Also in all cases the mappings 
were continuous. 

Inasmuch as a graph is a one-dimensional complex, we thus deal 
with continuous mappings of complexes onto each other. But in 
general a complex may be mapped onto another in many ways, all 
those many ways being divided into classes of homotopic mappings. 

The mapping which has been studied by us hitherto is suggested 
by biological considerations. We may call it, for lack of a better 
term, the natural mapping. We may ask the question: is the natural 
mapping the only possible for the complexes which we studied or 
is it, as we would expect in general, only a representative of a 
class of homotopic mappings. The answer to that question may be 
obtained only by a special investigation of all the transformations 
which have been studied hitherto. 

Consider the possibility of the existence of other mappings which 
are homotop to the natural mapping. In the natural mapping some 
points f;, fps fm» fpf, ate mapped onto the point f, of the primordial 
graph. Those points all represent different biological functions of 
the higher organism. They are the ‘‘partial biological functions”’ 
of the ‘‘biological function’’ of the primordial organism which is 
represented by the point f,. Other points fates ff of the trans- 
formed graph are mapped, under the natural mapping, onto another 
point fj of the primordial. If there exist other mappings, besides 
the natural, then some of the points of the transformed graph, say, 
f., 3, which are mapped onto fj under the natural mapping, may 
now be mapped onto /f,. 

What would be the biological interpretation of such another 
mapping if it is homotop with the natural? 

All the points of the transformed graph T which map on the same 
point f, of the primordial graph P represent different biological 
functions subsidiary to the general biological function represented 
by the point f,. Thus, for example, f, may stand for perception of 


light. Then f;,/,,/,, may stand respectively for perception of light 
of different colors. 
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Each biological function in a differentiated* higher organism is 
performed by a special type of cells best adapted for this performance. 

If now through a homotopic deformation of the mapping certain 
points, f, and f,, which under natural mapping correspond to fs 
now correspond to /,, then the natural biological interpretation of 
this would be the following: 

Through continuous biological changes, certain cell types, origi- 
nally performing biological functions subsidiary to one primordial 
biological function, now begin to perform biological functions of a 
different kind, corresponding to a different primordial biological 
function. 

It is a known fact that certain cell types which perform a certain 
biological function in one organism may begin to perform entirely 
different biological functions through gradual evolutionary changes. 
The character of those cells also undergoes a change. But ap- 
parently not any cell type may change into any other one, and the 
important problem arises as to the criterion which determines what 
kind of cell types can change into what others or, in other words, 
what kind of biological functions can be changed into other dif- 
ferent biological functions. 

Very tentatively we like to suggest the following topological 
criterion: 

If there exist two homotopic mappings of T onto P, such that in 
one of them the point f, of T maps onto the point f,, of P, and f, 
of T maps on f,, of P; while in the other either both f, and /, map 
on f,, or on f,,, or 7, maps on f,, and f, on f,,; then the cell types 
of the higher organism represented by 7, which perform the biologi- 
cal functions represented by f, and f,, can change one into the 
other. If no such homotopic mappings exist, no such change can 
occur. ; 

The criterion is in principle amenable to empirical verification, 
once the graphs T have been assigned to different higher organisms. 


This work was aided by a grant from the Dr. Wallace C. and Clara 
A. Abbott Memorial Fund of The University of Chicago. 
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CONTRIBUTIONS TO THE THEORY OF ACTIVE TRANSPORT: 
Il. THE GATE TYPE NON-CARRIER MECHANISM AND 
GENERALIZATIONS CONCERNING TRACER 
FLOW, EFFICIENCY, AND MEASUREMENT 
OF ENERGY EXPENDITURE 


CLIFFORD S. PATLAK 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO* 


The gate type non-carrier mechanism, an active transport model, is 
discussed. In this mechanism, the actively transported particle passes 
through the gate itself by means of a series of reorganizations of the -ac- 
tive transport mechanism. The net rate of transport, the rate of transport 
in either direction, and the efficiency of this model are analyzed. It is 
shown that on the basis of these analyses alone, this mechanism cannot 
be distinguished from a carrier mechanism, 

Three generalizations which apply to many individual type active 
transport models are then discussed. ‘These pertain to (1) the depend- 
ency of the flow in one direction on the concentration of the particles on 
the opposite side of the membrane, (2) the possibility of very high effi- 
ciencies for these models independent of the rate of the active transport, 
and (3) the methods whereby the energy expended in the active transport 
may be experimentally found. 


Introduction. In a previous paper (Patlak, 1956), a classifica- 
tion of the different types of active transport mechanisms was 
made. The possible mechanisms were first divided in two major 
groups. The first was called a general type of active transport 
mechanism. If A is the actively transported particle, then in this 
case the mechanism acts on all of the A particles simultaneously. 


*This research was supported by the United States Air Force through 
the Air Force Office of Scientific Research of the Air Research and De- 
velopment Command under contract No, AF 18(600)-1454. Reproduction 
in whole or in part is permitted for any purpose of the United States 


Government, 
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The second group was designated as an individual type of active 
transport mechanism. For this group, the active transport mecha- 
nism acted on one or a few of the A particles at a time. This may 
be accomplished in two general ways. The first is to carry the A 
across the active transport region via a carrier. The second way 
is to allow the A to cross the region by itself under more simple 
conditions. This will be mediated by a non-carrier mechanism. 
Although the general type of active transport mechanism does not 
use a carrier, we shall limit the name non-carrier mechanism to the 
individual type of non-carrier mechanisms. With a little care, this 
should cause no trouble. (In the previous paper, this mechanism 
was called a Maxwell Demon. This name was chosen because the 
method of operation of the Demon, as discussed by Maxwell, was 
an example of the non-carrier method. However, the name ‘‘Max- 
well Demon’”’ signifies not just a particular transport mechanism, 
but a mechanism which can violate the Second Law of Thermody- 
namics. As stated explicitly in the former paper, we did not en- 
visage any Violation of this Law and in fact required our Demon to 
have an energy input coupled with itself. However, because of 
this association of Maxwell’s Demon with the Second Law of Ther- 
modynamics, we shall no longer use this name.) 

The non-carrier mechanisms may themselves be further classi- 
fied using the methods of transport involved. Thus, they may act 
by: (1) increasing the velocity of the A particle so that it crosses 
the active transport region in a shorter time; (2) modifying the A 
particle so that it passes over the potential barrier more readily 
either by (a) increasing the energy of the A or (b) altering the form 
of the A so that the membrane becomes more permeable to it; and 
(3) decreasing the resistance of the region to the flow of the A by 
either acting on the region (a) immediately adjacent to or (b) at a 
distance from the non-carrier mechanism. Obviously, increasing the 
velocity of the A particle also increases its energy. However, 
category 1 is listed separately because its mode of action is con- 
ceptually different from that of just increasing the energy. Non- 
carrier mechanisms may also be combinations of more than one of 
the above categories. Some examples have been discussed in the 
previous paper. 

As discussed in that same paper, the non-carrier mechanisms, 
unlike many of the carrier mechanisms, do not operate over large 
distances, but are confined to perform their action over distances 
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of the order of a few Angstrom units. However, for an active trans- 
port mechanism to be effective, it need not act over the entire 
membrane but only across any potential barrier such that the per- 
meability of the barrier to passive diffusion is small. In this way, 
the short distance of action of the non-carrier mechanism does not 
disqualify it as a possible mechanism. The thickness of the mem- 
brane across which active transport takes place is usually of the 
order of 100A. Thus, insofar as A is concerned, we may conceive 
of this membrane as being composed of many pores through which 
A can diffuse. The non-carrier mechanism may be placed at some 
point within many of the pores and in many cases may be an inte- 
gral part of the membrane. Thus the potential barrier mentioned 
above is either the non-carrier mechanism itself or is an ‘‘obstruc- 
tion’’ in the pore close to the non-carrier. 

We shall now analyze one type of non-carrier mechanism. After 
that, some general conclusions which are applicable to many dif- 
ferent mechanisms of the individual type of active transport will 
be discussed. 


Gate type non-carrier mechanism. In order that a non-carrier 
mechanism may decrease the resistance of its immediate surround- 
ings to the passage of A, it must act like a ‘“‘gate.’’ Thus we will 
call this mechanism the gate type non-carrier mechanism. The 
gate reacts with A on one side of a potential barrier to form a com- 
plex which then, without any translational movement, undergoes a 
reorganization such that in its new form A can leave the complex 
on the other side of the potential barrier. The gate then undergoes 
a second reorganization so that it returns to its original form. The 
energy source for this series of reactions may enter into any of the 
steps. 

A particularly simple example of the above mechanism is repre- 
sented in Figure 1. In this case the initial form, G,, of the gate 
is represented as a molecule which has an ‘‘opening’’ on the left 
side. Some examples to illustrate what the term ‘‘opening’’ may 
represent are: (a) the type of molecule present in the theory of the 
‘“‘lock and key’’ type of reaction; (b) a long molecule to which the 
A may attach itself at some center point; and (c) a portion of a 
molecule with which the A forms a chelated complex. The A pres- 
ent on the left side, A), can react with G, to form the complex G,. 
(The terms ‘‘left’’ and ‘‘right’’ are used instead of the conventional 
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FIGURE 1, Scheme of a gate type non-carrier mechanism. See text. 


terms ‘‘inside’’ and ‘‘outside’’ in order to emphasize the fact that 
the A which reacts with the gate is the A which is immediately ad- 
jacent to it and not the A which is in the solutions on either side 
of the membrane.) The complex G, then undergoes a reaction in 
which the molecule ‘‘closes’’ on the left side and ‘‘opens’’ on the 
right to form G,. This product in turn then reacts to release the A 
to the right side, A, and leaves the gate in the form G, which can 
react again to form G@,, (It should be emphasized that G,, G,, G ss 
and G, are all at the same point within the membrane. That is, 
the arrows in Figure 1 denote chemical reactions and not physical 
displacements.) For simplicity of calculation we will assume that 
the kinetics of each step are of the simplest type, with the reaction 
rates as shown in Figure 1. The necessary energy for the active 
transport arises by the further reaction of substrates with any one 
or more of the four steps. The substrate, S,» reacts with G, (and 
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A, forg = 1) toyield Gg 1) and a product, P ,» land A, for(g + 1) = 4], 
as shown in Figure 1. The reactions are ahown in the given direc- 
tion because we are assuming that the active transport is from left 
to right. It must be emphasized at this point that we are making 
no claims whatsoever that this particular model is applicable to 
any biological case. However, it does exhibit many of the proper- 
ties which many carrier mechanisms possess and cannot be dis- 
tinguished from a carrier mechanism on the basis of the measure- 
ments of the flux of A alone. Thus, this model is being treated 
not only to illustrate certain features which will be discussed more 
generally in the next section, but also to emphasize that there are 
different active transport models besides those of the carrier type 
which may be used to explain experimental results. 

We shall adopt the customary convention of omitting partition 
and activity coefficients and speak only of concentrations. Further, 
the A and the other molecules involved in the active transport will 
be assumed to be electrically neutral. These assumptions will be 
discussed at length at the end of this section. 

Only the steady-state system will be discussed. The rate of the 
reaction, which is equal to the net rate per unit area that A is 
transported across the membrane via the active transport mecha- 
nism, will be denoted by J as in the previous paper. The concen- 
trations of the different components will be expressed by placing 
parentheses around their names. The concentration of the A, S’s, 
and P’s will be in mole//, while the concentration of the gate in 
its different forms will be in mole/cm?. By making the following 
definitions, 


w, =k, +k, (S,)] (A); 

w,=lk, +h, ,(S,)], (9 = 213,4)5 (1) 
73) 3 < [K_, An oe 3 (P,)] (A,); 
w_, = [k_, + (Ps yl, (g =a 1,2,4); 


ks 2 
the rate of reaction can be written as 


J=w,(G,) — w_,(@,); 
= W_(G,) - w_,(@5); (2) 
= w,(@,) - w_ (G4); 
= w,(G,) — w_,(G@;,). 
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Denoting by G, the total concentration of the gate in all of its 
forms, i.e., 


Gr = (G@,)+(G,) + (G,) +(G,), (3) 
we find from the above equations: 
Jem 
4 4 
a! Pit I w_, 
Gy - g=1 g=1 : 
W,W,W,W, W_,W_,W_,v_, 
Des a. jp oret mn + Wy 54. 0)(cg4.1) + M—ces1)) 
g=l rr 
(4) 
where W, = W,, W_, = W_ , and w_, = w_ 


We sBall now discuss equation (4). We first see that the choice 
as to which S,’s and P,,’s are zero or non-zero, although affecting 
the mechanism of the active transport profoundly, does not alter 
the form of the resultant rate equation. If there is at least one S 
and P, equation (4) becomes for large values of the A’s, 


J=K Z a 5 
(A) (4)|’ _ 
where K and kK’ are positive constants with Kk’ 4 1. 

If there were no active transport, i.e., all of the S’s and P’s are 
absent, K’ would equal one and equation (5) would reduce to the 
type E of T. Rosenberg and W. Wilbrandt (1955). It can be easily 
shown from equation (4) that dJ/d A, > 0 and dJ/dA, < 0; both re- 
sults are expected from an examination of the model. Therefore, 
we see that the graph of J vs. (A;) with constant (A,), and J vs. (A,) 
with constant (A,) would be as shown in Figures 2a and 2b. Also, 
for a constant ratio, (A,/A,), the graph of J vs. (A;) would be as 
shown in Figure 2c. The size of the units for both the J and the 
A’s depends upon the particular values of the w,.’s. The graphs 
of Figures 2a and 2b are to be expected from the model. However, 
the result of J approaching zero as A, and A_ both become large, 
as shown in Figure 2c and equation (5), may not be so obvious. 
The explanation of this result is that when the concentration of 
the A’s is large, for each G, that yields an A,, another A. instantly 
combines with the resultant G, to form a G,. An analogous ex- 
planation applies to the A), with G, and G, in place of G, and G,. 
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FIGURE 2a. Graphs of equation (4), The X’s, which may be easily de- 
rived by suitable rearrangements of equation (4), are functions of the 
Wig» &T, (A,), and (Aj). Figure 2a is for J vs. (A) with (A,) held con- 
stant; 2b is for J vs. (A,) with (A;) held constant; and 2c is for J vs. 
(A) with (Aj)/(A,) held constant. The figures are drawn for the case of 
J>0. 


Thus, all of the gate is in the form of G, and G,, and therefore there 
is no net transfer of A from left to right. 

We may now use our model of the gate type non-carrier mecha- 
nism as a particular example in a discussion of the flow of tracers. 
Using an analogous notation to that used in the previous paper, we 
represent by J and J the flow of A per unit area and time across 
the active transport mechanism from the left to the right side and 
the flow from the right to the left side respectively. Further, let 
—» or « above the symbol of a molecule indicate the molecule 
which is combined with the tracer of A, or A, respectively. Thus, 
we have 


j= w(G,)- w_,(G,) = w. (Go) = w_, (Gs) =W, (2); 

5 = w_, (Gq) = wg (G,) - wy (G@,) = v5 (G,) ~ 4 (Gy); 
Jud —J=0,(G,)—w_,(6,); (6) 
(4) = (@,) + (@,); 

(G4) = (@,) + (G5) 
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FIGURE 2b. See legend for Figure 2a, 
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FIGURE 2c. See legend for Figure 2a, 
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p) 


G 
hal’ B 
ihe D W WeW, (wW_,wW_,W_, + W_,W_,w, + W_,W,W, + WyW,wW,) 
(8) 
G 


W_,W_,W_,(w, Ww, + WeW,W_, + WyW_,wW_, + W_,W_,W_,). 

The graph of 7 vs. (A,) for constant (A,) is shown in Figure 8a. 
Since w, is a constant -(A;) and w_, is a constant -(A.), we see 
that J and J are functions of both A, and A, unlike the case of 
Fick diffusion. However, depending upon the values of the reac- 
tion rate constants, 0J/d(A,) may be approximately zero and thus 
7 may for all practical purposes be independent of A,. An analo- 
gous situation will also apply to the dependence of J on A,. For 
example, if w, >> w_, and w, and w, are greater than w, and w,, 
we have the following three cases. In the first case, if w, >> w, 
>> w_,, then ‘J = constant - (A,) and is independent of A;. In the 
second case, if w, >> w, >> w_, (i.e., the reaction A, + “4 + G, 


is the rate limiting reaction in the left to right direction), J = con- 
Ky Kins 
K3=1 
a 
K'=10 
fe) 1 


1 re 


FIGURE 3a. Graphs of equation (8). The K’s, which may be easily de- 
rived by suitable rearrangements of eguation (8), are functions of w4,, 
Gr, (A,), and (A). Figure 3a is for J vs. (Aj) with (A,) held constant; 
and 3b is for J vs. (A,) with (Aj) held constant, The restriction that J > 0 
is not used for these graphs. 
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stant -(A_) (A,) and thus 7 is directly proportional to A). In the 
third LA8e; if wy << w, << w_,, then J ~ constant: (A VA) and 
thus J is Hreeaete S oporicne! to A; A similar type of depend- 
ency of Jon (A) for constant (A)) set be seen in Figure 3b. 

We shall now turn our attention to a consideration of the effi- 
ciency of this mechanism. By the term ‘‘efficiency’’ we shall mean 
the usual thermodynamic definition: the ratio of the energy gained 
by the system to the energy put into the system. Thus, the effi- 
ciency of the model will be independent of whether the energy put 
into the mechanism is a byproduct of a metabolic reaction which 
performs primarily some other function besides the active transport, 
or whether the energy is from a reaction the main function of which 
is to supply the energy for the active transport. It will also be in- 
dependent of the amount of A which crosses the membrane without 
passing through the active transport mechanism. This follows from 
the fact that we are looking at the efficiency of the active trans- 
port mechanism and not the efficiency of the entire system. 

For this calculation, we must consider what happens to the sys- 
tem as it passes through a complete ‘‘cycle.’’ Thus, for each 
particle of A which is actively transported from the left to the 
right side by any means whatsoever, we must also consider the 
number of particles of A which are transported, with or without an 


> 


O° kK," 
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FIGURE 3b, See legend for Figure 3a. 
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expenditure of energy, in the same or opposite direction by all 
other means involving the active transport mechanism. For each 
of the pathways, a certain amount of energy (positive or negative) 
will be given to the system by the transfer of the A’s from one side 
to the other, and a certain amount of energy will be expended by 
the system to accomplish the transport. As stated before, the 
ratio of the energy given to the energy expended, for all of the 
pathways, is the efficiency of the mechanism. 

For the particular mechanism under consideration (Figure 1) we 
see that A may be transferred from one side to the other either by 
means of energetic pathways (involving one or more of the S’s and 
P’s) or by a non-energetic pathway (involving no S or P). The frac- 
tion of the particles of A which are transported by any particular 
pathway may be readily calculated. However, the result of this 
calculation is extremely complex. The presence of an energetic 
reaction usually implies that the number of particles of A trans- 
ported by the parallel non-energetic pathway is negligible com- 
pared with the number transported by the energetic pathway itself. 
Thus, due to both of these reasons, we shall make the simplifying 
assumption that, where there is an energetic reaction, the parallel 
non-energetic pathway may be ignored. That is, where kegs #0, 
we shall assume that both k, and k_, are negligible compared with 
ky oe Thus, the energy put into the mechanism would be the sum 
of the energies dissipated by the energetic reactions. The overall 
reaction of the model may then be written as 


Agta + Dodi bg + 8,:¢2 Art Py + Pot Py t Pa, (9) 


and the energy introduced into the system would be the sum of the 
differences of the partial molar free energies of the P’s and the 
S’s. If we let u, represent the standard chemical potential, then 
the energy introduced into the system under the assumptions made 
earlier in this section is given by 


: RT lo (S,) + -p ; (10) 
pe ° (P,) o,S, 0,P, 
g= 


The energy derived from the mechanism is simply the energy 
gained by the system when one mole of A is transported from the 
left to the right side which is the difference between the free 
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energy of A_ and the free energy of A}. Thus, under the assump- 
tions listed before, the energy gained is given by 
(A,) 
(A;) 


Hence, the efficiency, which is the ratio of expressions (11) and ; 
(10), is 


RT log (11) 


(A,) 
A 
efficiency = idea Ts eee (12) 
Dn RT 108 (By P ) ~ Ho,S, - mA 
g=1 & 


The principal function of an active transport mechanism is to 
transport A from a lower to a higher chemical potential. As dis- 
cussed in the previous paper, this need not always happen, but for 
simplicity we shall assume this and thus 


(A,) > (A)). (13) 


If we did not make this assumption, the efficiency, as defined 
above, would be negative because no energy would be gained by 
the transport of A from left to right. For that case, the purpose of 
the active transport mechanism would be to transfer A from left to 
right more readily than otherwise and not to accumulate A on the 
right-hand side. The free energy change of the reaction from left 
to right of equation (9) must be less than zero because otherwise 
the reaction would go in the opposite direction which would mean 
that the mechanism would be pumping A from right to left. Under 
the assumption of equation (13), the efficiency in equation (12) is 
positive and less than or equal to 100%. 

To transform equation (12) into a more workable form, we intro- 
duce the subscripts eg, to denote the equilibrium values and make 
use of the fact that 


(od) aon, AO eae 
D) (Ho,s, ~ Hop) = RT log ett + RT Donde geet (14) 


leq eal 
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By equilibrium we mean equilibrium with respect to the reaction 
shown in equation (9), which may be found by letting J = 0 in equa- 
tion (4), where kis. #0, with k= kh, =(. Thus, 


(A,) 
(A) 


aes) at es T(S\(P 22 
l 2€Q & &,e€ 
eocamey PS Bees ge 


log 


efficiency = (15) 


need 


We see that when the system approaches equilibrium the efficiency 
of the pump approaches 100%. 

As seen from equation (15), if (A)) is increased or (A_) is de- 
creased, the efficiency will decrease. However, as discussed pre- 
viously, the above changes also cause an increase of J. There- 
fore, increasing the value of J for a constant A_or A, causes a de- 
crease of the efficiency. Thus, a 100% efficiency holds only when 
J=0. For a finite value of J, the efficiency, as expected, is less 
than 100%. As seen in Figure 2c, by keeping (A,/A,) constant we 
may have an increase of J with a decrease of A). Thus, it is pos- 
sible for an increased J to have an increased efficiency. 

The actual relationship between J and the efficiency may be 
found from equations (4) and (15). However, it is not only extremely 
difficult to solve the two equations simultaneously, but the result 
does not appear to be of any particular value. What is of most in- 
terest is that for a given value of the velocity of the reaction we 
can get any desired value for the efficiency, between 0 and 100%, 
by a suitable choice of parameters. Contrawise, for any given ef- 
ficiency, we can have any value of J. Thus, if in an experimental 
situation the efficiency is fairly high (70 to 80%) and the velocity 
of an apparent active transport is also high, there is no apriori 
reason to assume that there is no active transport. Although it is 
very improbable that such high efficiencies will be observed in bio- 
logicai material, especially in conjunction with high rates of trans- 
port, it is not impossible in theory. 

To conclude this section we shall briefly discuss the assump- 
tions made earlier that the partition coefficient is the same on the 
left and the right side, that the activity coefficient is equal to one, 
and that all particles are electrically neutral. 
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In many discussions of active transport in the literature, it is 
envisioned that A is taken from the solution on one side of the 
membrane and transported to the solution on the other side. How- 
ever, as pointed out earlier, this is not even theoretically possible 
for non-carrier mechanism. In this case, the A, and the A, are both 
within the membrane, and, unlike the former examples where the 
partition coefficients of the A in the two solutions are usually very 
similar, it is less probable that the partition coefficient for A, and 
A, will be the same. Thus, the omission of the partition coeffi- 
cient is an assumption, the validity of which is open to question in 
any particular case. However, the non-validity of this assumption 
does not affect the conclusions of the discussion on the efficiency 
of the active transport mechanism. 

In many experimental situations a charged particle is actively 
transported. If A is a charged particle, the assumption that the ac- 
tivity coefficient equals one is of dubious validity. However, due to 
the presence of other electrolytes, the activity coefficient of the A 
may be independent of the concentration of the A. In this case, 
the partition and activity coefficients may be included in the reac- 
tion rate constants and thus the form of equations (4) and (8) will 
not be altered. Since the values of the reaction rate constants are 
very rarely known, the inclusion of the partition and activity coef- 
ficients in the reaction rate constants will therefore cause no 
trouble. 

The introduction of charged particles complicates the problem. 
We have here a situation in which the components of the reaction 
do not all remain in the same environment, unlike most of the chemi- 
cal reactions which are considered in the literature. Instead, as 
seen from Figure 1, part of the reaction occurs on the left side, 
part in the region between the two sides, and part on the right 
side. Therefore, the introduction of an electric potential and 
charged components will affect the various reaction rates. If there 
is an electric field, not only must there be very specific assump- 
tions about the charge and position of the various components, but 
the position of the activation energy hump must also be made ex- 
plicit in order to determine the values of the reaction rate con- 
stants (Glasstone et. al., 1941; Karreman, 1951). Also, since the 
non-carrier mechanism is of the order of only a few A in width, 
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then, as discussed in the previous paper, electric neutrality can- 
not be assumed. 


Generalizations. In the previous section, certain results were 
found for the gate type non-carrier mechanism. However, some of 
these results are valid not only for that model, but also for many 
of the active transport mechanisms that have been proposed in the 
literature. Three generalizations of some of those results will be 
discussed below. 

1. The first generalization is concerned with the dependence of 
the flow of tracers on the concentration of the A’s on the ‘‘oppo- 
site’’ side. That is, if j represents the net flow of A across the 
membrane per unit time and area and the subscript s refers to the 
solutions on either side of the membrane, then 7 may depend on A,. 
and ) may depend on A, 

We shall show how this may be so under very general considera- 
tions without the necessity of assuming a change in either the mem- 
brane permeability or any chemical reaction. Let us consider the 
situation where A is transported by an individual type of active 
transport mechanism. The A will be actively transported by the 
aid of one or more ¢ransporting-particles. In this mechanism, the 
A; through a series of chemical reactions, combines with a trans- 
porting-particle in its ‘‘initial state’’ on the left side of the active 
transport mechanism. (In the discussion that follows, we shall 
focus our attention on this first transporting-particle in order to 
explain the possible dependency on the (A_) by the rate, J, at 
which A, is transported from left to right), This combination will 
undergo further chemical reactions and physical displacements (as 
in a carrier) that culminate in the A, being released from the trans- 
porting-particle which will then be in its ‘‘final state.”’ (In a sim- 
ple carrier model the final state of the transporting~particle is the 
carrier on the right side of the active transport mechanism; in the 
gate type non-carrier model discussed in this paper, the final state 
of the transporting-particle would be G,.) The transporting-particle 
in the final state will return to its Fnitial state either spontaneously 
or through a second series of chemical reactions. 

If we assume that A is transported from left to right and from 
right to left by the same transporting-particles, A will be trans- 
ported from right to left by reactions which are the reverse of the 
above. Many specific examples of the above general mechanism 
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may be found in Rosenberg and Wilbrandt (loc. cit.). The net ve- 
locity of the reactions will be in the first direction due to the 
energetics of the chemical reactions. For a given concentration of 
substrates, products, and A,, and for a given V, and assuming that 
there are no further transporting-particles, J will be a monotoni- 
cally increasing function of the concentration of the transporting- 
particle available to react with the A,, i.e., the concentration of 
the transporting-particle in its initial state on the left side of the 
transporting region. However, this concentration will in turn be a 
function of the rate at which A, is transported from right to left, 
x This is because, by -deatite with A, the transporting-particle 
usually goes from its final state to its initial state at a different 
rate than the transporting-particle would go by itself. Thus, if the 
rate with the A, i.e., the ‘‘back’’ reaction, is larger than without 
the A, then in the steady state with the presence of A, there will 
be a greater concentration of the transporting-particle in its initial 
state than if there were no A. An extreme example is when the 
transporting-particle is normally destroyed in its final state and 
thus its rate of return to its initial state is zero. Since Jisa 
monotonically increasing function of (A,), then J will be a mono- 
tonically increasing function of (A,). Contrawise, if the rate of 
transformation of the transporting-particle from its final to its ini- 
tial state after reacting with A, is smaller than its rate without the 
A,, then by a similar argument as befgre, J would be a monotoni- 
cally decreasing function of (A,). If J is independent of (A_) be- 
cause some other step in the transport process is the rate limiting 
reaction, then 7 will also be independent of (A,). Examples of 
these three cases have been discussed in the paragraph following 
equation (8) for the dependancy of J on (A)). 

The resultant J could have a complex dependency on the A, which 
is not a simple monotonically increasing or decreasing function of 
the (A,), This is so for at least two reasons: First, if there are 
further transporting-particles involved in the active transport, an 
analogous argument would apply to each of them. The second 
reason is that the presence of the A, may affect the rate at which 
the combination of A, and the transporting-particles undergo their 
further reactions, for example, when the transporting-particle re- 
acts with more than one A at a time. (The general conclusion for 
J can also be found using other mechanisms, such as that pro- 
posed by L. Kirschner, 1955.) 
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By a completely analogous discussion, we may arrive at the 
same conclusions as above with regard to the dependency of J on 
(A)). Since ; and I, 7 and Fi A, and A,, and A, and A are re- 
lated, j may be a function of A_ Oe well as of A), etd x may be a 
function of A, as well as A re" Hor the case of epuuess with the 
same Anginene discussed above, but with the energetics of the 
reactions such that there is no active transport, as in ““exchange- 
diffusion,’’ the above argument is also valid. If there were an ex- 
tremely large number of transporting-particles so that at any time 
the number available to react with the A _is independent of the con- 
centration of the A, then obviously the 7 would be independent of 
Zo. and j j would be independent of A). 

‘Tf the concept of the “Serer rae is further generalized 
to include a channel through which A passes, the mechanism of 
A. L. Hodgkin and R. D. Keynes (1955) is then included. (This 
discussion does not, of course, cover all aspects of their model, 
but only deals with the dependence of the j on A.) In the case 
of the usual type of passive diffusion (following Fick’s Law), 
since the molecules of A). and a can avoid each other as they 
pass across the membrane, the channels through which the mole- 
cules pass—which are of course hypothetical ones and part of the 
actual channel—are effectively infinite in number. If the concen- 
tration of the A’s becomes extremely large so that the A, and A, 
would collide with each other as they passed through the feria 
we would then have, in effect, a finite number of channels and the 
situation would reduce to that discussed by Hodgkins and Keynes. 

From the analysis of the non-carrier model of the previous sec- 
tion and from the discussion above, some further facts about j and 
j may be noted. It has been customary to assume that if A, =A, 
and if, 7z ji, then in most cases j represents passive He 
while 4 —j represents the active transport. The only qualification 
to this seems to be that of H. H. Ussing (1949) who stated that 
some of the j might be due to ‘‘exchange-diffusion.’’ However, as 
we have seen in our model of the gate type non-carrier mechanism, 
and as is obviously true for any type of active transport pump, 
there is some flow in both directions due to the action of the pump. 
The amount being transported in the ‘‘back’’ direction (in this 
case J) would depend on the concentration of the A’s as well as 
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on the parameters of the pump and may be negligible for many ac- 
tual cases. However, the value of J approaches that of i; when 
the values of A, and A, approach the equilibrium values for the 
pump. Thus, the presence of a 7 not equal to zero may indicate a 
backward passive diffusion, it may be due to an ‘‘exchange-diffu- 
sion,’’ or it may simply arise as part of the operation of the active 
transport pump, or any combination of the three. In any case, this 
must be considered a possible cause of j j and should be tested for. 

9. The second generalization is concerned with the possibility 
of 100% efficiency for the active transport mechanism. Before dis- 
cussing this further, we shall discuss some necessary assumptions. 

We shall consider active transport mechanisms of the type dis- 
cussed in the previous generalization and shall first assume that 
there is only a single pathway for the A’s to traverse the active 
transport mechanism. The overall chemical reaction of the system 
is then 


mS, + mS, +ere+ mS, + EA, =? 2,P,+ 
MoPyte+++N Pi + €A,,°°(16) 


where the m,, n,, and € are constant coefficients and where the g 
need not be 4, as in the previous gate model, but may have any 
value depending on the particular model. For the case of the model 
discussed in Rosenberg and Wilbrandt, (Joc. cit.), g=2. We may 
note that equation (16) will be valid for many types of carrier and 
non-carrier models. If the rate of the active transport mechanism is 
zero, then due to the above assumption, all of the constituents of 
the mechanism will be in equilibrium. 

We will further assume that A has a valence, 2, and that there is 
an electric potential, V (equal to V_ - V,), across the active trans- 
port mechanism. The next necessary assumption is that either all 
of the S’s and P’s are electrically neutral or, if not, then those 
that have an electric charge do not move during the transporting 
reaction. Finally, we shall make the simplifying assumptions that 
the activity and partition coefficients for the A, and A, are the 
same and that the activity coefficients of the S’s and P’s remain 
constant as the concentration of the S’s and P’s change so that the 
activity coefficients may be included in the »,’s. As in the previ- 
ous section, the validity of the possibility of 100% efficiency is 
not affected by the non-validity of the latter assumptions. 
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The efficiency of the active transport mechanism is defined in 
the same manner as before. By an analogous argument to that 
used in the previous section, we have 


(4,) 
(A)) 


ig f Sz)he 
& 
ei, Pr tn, S (Pyne j* Metos, "Par, 


As stated above, when the rate of the active transport mechanism 
is zero, all of the components of the mechanism will be in equilib- 
rium. Thus, the reaction of equation (16) can reach an equilibrium 
state. (This is contrary to the Franck=-Mayer model of a solute 
diffusion pump (1947) where, when the rate of the active transport 
is zero, the carrier mechanism is not in equilibrium and then the 
analogue of equation (16) for this case would not go to equilibrium. 
This may be explained as due to the violation by the Franck=Mayer 
model of the first assumption above. This is in agreement with S. 
Golden (1948) who implicitly stated in his correction to the Franck- 
Mayer paper that 100% efficiency can only be achieved when this 
assumption is valid for their model). The same type of argument 
used in deriving equation (14) can now be used, and equation (17) 
becomes 


RT log - + 2FV 


efficiency = 
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As was mentioned in the previous section, equation (16) differs 
from the usual chemical reaction equation in that all of the reac- 
tion does not take place in the same region. It is for this reason 
that the 2FV term appears in equations (17) and (18). From the 
equation above, we see that as the system approaches equilibrium 
the efficiency approaches 100%. Hence, the same argument as be- 
fore, regarding the possibility of high efficiencies for finite J’s, 
applies. 


(18) 


B,eq 
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It should be emphasized that in the equations (17) and (18) of 
the efficiency we have considered the A and the electric potential 
across the active transport mechanism—not across the entire mem- 
brane. The reason that the latter is not used is that if the system 
as a whole is not in equilibrium, then j 4 0 and the A diffuses from 
the solution to the pump on the left-hand side and from the pump to 
the solution on the right-hand side. Thus, there is a loss of free 
energy resulting from this diffusive flow. The energy gained by 
the pump’s action is shown in the numerator of equations (17) and 
(18). Hence, a smaller value will obtain if we use the values for 
the concentration of A in the solutions and the value for the elec- 
tric potential across the membrane, in the numerator of equations 
(17) and (18), instead of both of these values in the regions adja- 
cent to the active transport mechanism. This may also be seen if 
we note that, for 7 to be positive, (4; ) exp{(2F/RT)V, }must be 


greater than (A;)expi(2F/RT)V,}; and (A,)expi(zF/RT)V,}must 
be greater than (4, ) exp{(2h/RT)V, $. Thus, letting V’ be the 
electric potential between the left- and right-hand solution (equal 
to ate - Vie) (A, )/(A,,)I exp{(zF/RT)V’}will be less than 
[(A)/(A,)] exp {(2F/RT)V}. 

If there is passive diffusion across the membrane, then j 4 J. 
Thus, the pump cannot be in equilibrium when the system is in 
pseudo-equilibrium, i.e., whenj = 0,J #0. Therefore, (A)) #(A, e8 
and (A) # (A, .,) even though the A, and A, will be in equilibrium 
with the A). and A, respectively. 

Instead of measuring the efficiency of the active transport mech- 
anism, we may wish to find the efficiency of the system as a whole. 
Under the assumptions listed above, we see that for each molecule 
of A which is pumped across the active transport mechanism, j/J 
will traverse the membrane. Thus, the efficiency of the total sys- 
tem is given by 


efficiency = 
(4,,) 
RT log —*~+ 2FV’ 
cs (Ai,) ei 
A m : 
RT log Seta 2FV, , pial = log (S,) eh cy 
(ar eal Ce (PF ae 
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3. The third generalization of this section concerns the measure- 
ment of the work expended by the system to accomplish the active 
transport. We shall first discuss the equation used by K. Zerhan 
(1956, p. 311) in order to achieve this measurement. Although he 
does not give the derivation of this equation, it follows directly 
from equation 1.2.7 of H. Linderholm (1952), Letting 2=+1, and 
¢ be the ‘‘active transport potential,’’ the equation of Linderholm 
is, in our notation, 


> 
j (Al) -#RW'+h 
; (A,,) 
Thus, the energy expended by the active transport mechanism, W, is 
=> 
j A 
W=-F6=RT loge + ot @ rs rs). py, (21) 


j (A;,) 


all terms of which are capable of measurement in many experimen- 
tal situations. However, as shown in the previous paper, the equa- 
tion from which equation (20) was derived is not applicable to car- 
rier or non-carrier models because the active transport potential 
applies only to the general type of active transport. For this latter 
type, the extra flow at any point due to the action of the active 
transport potential is equal to -uF (A)d¢/dz, where u is the mo- 
bility of A. However, for the individual type of active transport, 
the flow due to the action of the pump is equal to J , a constant 
which is usually dependent on the concentration of the A only at 
the two boundaries of the pump. Of course, for any particular ex- 
periment, a suitable value of ¢ may be chosen so as to yield the 
correct results. However, this value of ¢ will be the value for a 
‘mythical’? pump which performs its action via a general type of 
active transport mechanism to yield the required results. Thus, 
the value of ¢ will yield the energy of this mythical pump and not 
the value of the energy of the actual pump. (As a corollary of this 
argument, we should not expect the value of ¢ to be constant as 
the conditions of the experiment are changed by altering concen- 
trations and potentials, even though the pump itself is not basi- 
cally altered). Since most of the active transport mechanisms that 
have been formulated and applied to specific cases are of the in- 
dividual types, it would appear that equation (21) is of extremely 
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limited applicability. However, as will be shown, equation (21) 
can be used in many of these cases under certain conditions. 

The energy expended by an active transport mechanism of the 
type discussed under the second generalization is shown in the de- 
nominators of the right side of equations (17) and (18). The major 
difficulty in ascertaining the value of this energy is not only the 
difficulty of measuring the values of the S,, P,, and y,’s, but in 
even determining what the actual reactants are. Therefore, we 
would like to transform the equation into a form whereby some 
other variables which are measurable can be introduced instead of 
the above variables. To do this, it is necessary to define our sys- 
tem precisely. 

We will assume that the system is the same as that discussed 
under the second generalization, with the addition of two further 
assumptions. The first is that during the course of any experiment 
on the system, the S, and Be remain constant in value. The second 
is that there is no passive diffusion of A across the membrane, 
i.e., A only crosses the membrane via the active transport mecha- 
nism, and therefore j= J. These two assumptions imply that the 
equilibrium values of the S’s and P’s are AS = oF and gh = a 
and if there is no net movement of A across the membrane, the 
chemical reaction given by equation (16) is in equilibrium. Also, 
at equilibrium (j = J = 0), A, and A, are in equilibrium with A, and 
A, respectively. Thus, the denominator of equation (18) may be 
written as 


W = RT log (Area) + aFV,, = RT log eee +2FV.,. (22) 
l,eq l,,eqg 

Hence, if an equilibrium value of A,, and A,, could be found and 
the electric potential measured at that time, the energy expended 
can be found. This will be further discussed later in this section. 
However, it may prove difficult experimentally to find an equilib- 
rium value of A, and Aj. so that an alternative approach must be 
used. 


» 
Let us assume that we have measured j /j and found that over a 
region of values of A, and A, 
s 


(4,,)]* 


= constant. 
(A,,) P 
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where y is a constant. If we make the necessary assumption that 
there is at least one equilibrium value for A). and A_ within the 
range of values of A for which equation (23) is valid, then for this 
value, j/j=1. Therefore, the constant of equation (23) equals 
LAS oes WA yi 1 If we further assume that the electric poten- 
tial remains constant during the course of the experiment so that 
je F al ae we see that equation (22) may be written as 


(4,,) 
(4, ) 


pre 7 
= - loge + RT log + 2FV’, (24) 
For the case that the y= 1, the denominator reduces to equation 
(21). 

The assumptions which were needed in order to arrive at the 
conclusions of the last two generalizations are not very stringent 
and may be applicable to many actual cases. Most of them are 
straightforward and their validity depends strictly upon the partic- 
ular experimental situation. Hence we will not discuss them any 
further except for the last two assumptions: equation (23) and 
Vv’ =V*_. It is possible that these two assumptions will be valid 
for an arbitrarily given region of concentrations, although it is rather 
improbable in general. However, if we restrict the experiment to 
conditions of low concentrations of A). and A, then it is prob- 
able that, if the other assumptions are valid, not only will the two 
assumptions be valid but equation (23) will also be valid for at 
least one equilibrium value of A... and A, The constancy of V’ 
is due to the presence of other electrolytes. If the concentrations 
of the A’s are small, then the electric potential will be essentially 
determined by the other electrolytes and V’ will remain constant 
as the concentrationof the A’s are changed. (However, if the mem- 
brane behaves as an A-electrode, then an alternative approach— 
keeping the V’ constant by means of an externally applied potential 
—might be used). 

The validity of equation (23) arises from the fact that, at low 
concentrations of the A’s, the A’s will be the rate limiting factors 
in all of the reactions involving them throughout the membrane. 
Thus the j j and the j will be proportional to (A, Sd and (A, a re- 
spectively, where % will be the order of the rate limiting reaction 
in the membrane involving the A’s. Since V is constant, their ratio 
will satisfy equation (23). There obviously exists equilibrium 
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values of A, and A, which are small in value since from equation 


(16) it is the ratio of (A,)/(A,) which is important for equilibrium 
and not the actual values. Thus the necessary condition that equa- 
tion (23) be satisfied for equilibrium values of A, and A, will also 
be true. Therefore, in order to use equation (2th, we hoard limit 
ourselves to low concentrations, but, if this is not possible, equa- 
tion (22) may then be used. 

H. Ussing (1954, p. 419), discussed three independent experi- 
mental methods which could be used to find the value of ¢, and 
thus W if equation (21) was valid. These methods were based on 
the assumption that the active transport mechanism behaves as a 
battery which is in series with the frictional resistance of the 
membrane to A, both in parallel with the frictional resistance of 
the membrane to all of the other particles of the system (see his 
figure 3). Except for the replacement of the pump by a battery, 
this model is a representation of the assumption j = J which we 
have been using in this section. His first experimental method, in 
our terminology, is tolet A, = Ais and to choose V’ so that 7 = 0 
(i.e., V’ = V7), and to éqittts this V’ with ¢. [It is advisable to 
choose the A; and A_ sg that j is a maximum (see, for example, 
Figure 2c) sad thus the value of V’.. can be more accurately de- 
termined.] Ussing’s second ererimental method is to short- 
circuit the membrane (i.e., V* = 0), and to adjust A). and ae until 
j=0 [i.e., (A, /AI as (A WAL, eg dls and to equate ses log 


rs,eq 


[(A, MA, my with ¢. (These two experimental | methods are also 
suggested ‘by equation (21) since for j = 0, log G/5)-= 0.) From an 
examination of equation (22), to which equation (21) reduces for 
these cases, we see that both of the above experimental methods 
will yield the correct energy expenditure. Ussing’s third experi- 
mental method is to find @ via equation (21) with 740. As we 
have shown, this method will yield the correct energy expenditure 
only under special conditions. 

In an experiment on frog skin (H. Ussing and K. Zerahn, 1951), 
the first method led to an average value for & of 110 mV, which 
yields a W of 2550 cal/mole. The third method, on the basis of 
the data of tables 2 and 3 of their article, averaged 75 mV for dy 
which yields a W of 1750 cal/mole. Thus, the results of the two 
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methods do not agree, as might well be expected. [In another set 
of experiments on frog skin (L. Kirschner, personal communication, 
1957), the second method yielded a 4 of 157 mV while the third 
method gave a ¢ of 56 mV.] 

Equation (22) may also be used to test the validity of the as- 
sumptions that j= J and that there is no passive diffusion of A 
through the active transport mechanism. If these assumptions, as 
well as the others made in this section, are valid and if the con- 
centration of the S’s and P’s are not altered in the experiment, then 
W remains constant while the values of a Bee aS: Veg and 
the concentrations of the other particles in the solutions are varied 
separately. That is, the product 


(A, ae) FrVeg 
(ae) a2, 


remains constant for a membrane satisfying the above assumptions. 
This constancy is true even if the pump operates without trans- 
porting any A. That is, the active transport mechanism may work 
constantly. If an A particle is present to combine with the pump, 
it will be transported across. If the A is not present, the pump 
will continue to perform its activities as if the A were there and 
hence it will consume energy. For example, consider a slight 
modification of the gate model discussed in this paper. If the re- 
action S, + @G, =P, + G, also takes place, where now G, and G, 
are assumed not to conta an A, and the remainder of the reaction 
proceeds passively, we see that this_model will operate with or 
without transporting an A. 

Equation (16) will still be valid in describing the transport of 
the A and W will still be the work done by the active transport 
mechanism to transport one mole of A. But W will not be the total 
work done by the active transport mechanism during the time that 
the mole of A is transported. Thus, even for these cases, W will 
be constant and the product of expression (25) will remain con- 
stant. However, if there is a pathway for passive diffusion of A 
across the membrane, either parallel to the active transport mech- 
anism (j 4 J) or through the active transport mechanism, then the 
derivation of equation (22) is not valid. Hence, if the product of 
expression (25) is experimentally found to be non-constant as the vari- 
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ables of the system are separately altered, ‘‘eg’? now meaning the 
value of the variable when j = 0, this might indicate a passive diffu- 
sion pathway. Note that this result is contrary to equation (21), which 
was derived from the active transport potential concept without the re- 
quirement that there be no passive diffusion pathway. However, if 
the product of expression (25) appears constant, this does not in- 
dicate that there is no passive diffusion pathway across the entire 
membrane since the non-constancy of this expression might become 
experimentally obvious only for values of the variables in the non- 
physiological regions and therefore would be undetected. Also, 
since the equations governing the change of the product of expres- 
sion (25) are not known for the case of a passive diffusion path- 
way, no conclusion can be drawn from an apparent constancy of 
this product. 


The author wishes to express his gratitude to Mr. Peter Curran 
and Drs. Richard Durbin and Robert Macey for reading and dis- 
cussing this paper. 
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10 Rubles. 


Every mathematical biologist should read this interesting and stimu- 
lating book. Although the title seems to emphasize the experimental 
aspect of the problem, yet the book is thoroughly permeated with the 
spirit of mathematical biology and stresses everywhere the need for a 
parallel development of theory and experiment. True enough, the mathe- 
matics is rather elementary, and the reader will not find any such elabo- 
rate theoretical developments as are found, for example, in the pages of 
this journal. But wherever the author deals even with a very elementary 
mathematical treatment of a problem, the mathematical reader finds sug- 
gestions and inspirations for a more thorough analysis, 

Perhaps this is the appropriate opportunity for the reviewer to state 
publicly an opinion which he has held for many years. The use of com- 
plicated mathematics is not necessarily always a sign of a truly mathe- 
matical spirit of the approach. An example taken from a different field 
may illustrate the point. A pianist may possess a phenomenal technique 
but be relatively devoid of feeling in his performance. Another may show 
much more feeling, but have a poorer technique. The ideal musician 
must, of course, possess both feeling and technique. Similarly, a mathe- 
matical biologist, or for that matter any applied mathematician, must 
strive for both the use of the most powerful mathematical methods and 
for the understanding of the real spirit of a mathematical approach to 
natural sciences, Without this spirit his work may reduce to a mere ex- 
ercise in solving highly complicated and difficult equations, 

Though Professor Ivlev does not explicitly exhibit any complicated 
mathematical technique (which, of course, does not mean that he does 
not possess it), yet he clearly shows the understanding of the spirit of 
mathematical biology. 

We shall give only a couple of examples of the type of problems and 
their treatment as found in the book. 

After two introductory chapters, the author discusses in Chapter Ill 
the effects of concentration of food on the rate of its consumption, In 
the classical studies of A. J. Lotka and of V. Volterra, the total rate of 
consumption of food is considered as proportional to the product of the 
number of consuming individuals and the concentration or amount of avail- 
able food. This is equivalent to assuming that each individual consumes 
food at a rate proportional to the concentration or amount of the latter. 
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The author quite correctly points out that such an assumption may be re- 
garded at best as a first approximation for very low concentrations. Ap- 
plied to the whole range of possible concentrations it leads to a biologi- 
cal absurdity, namely, to the possibility of an animal consuming at any 
rate, no matter how large. The rate of food consumption is actually 
limited by a number of physiological factors, such as speed of swallow- 
ing, duration of stay of food in the stomach, etc., factors which may 
have nothing to do with the outside concentration of food. 

The author argues that with increasing concentration p of food equal 
increases of concentration should produce an increase of the rate r of 
consumption in proportion to the difference R ~ r between the physiologi- 
cally highest possible rate of consumption # and the actual consumption 
r, In other words, denoting by € a constant, the author postulates the 
equation 


dr 
ae E(R Ys (1) 
dp 


which upon integration gives 
r=R(1- eS P), (2) 


The author’s own experiments confirm equation (2) to a remarkable de- 
gree. From the experimental r,p curves the author calculates the values 
of € Those values vary from case to case, but are all of the same order 
of magnitude, 

The fact that the simple relation (2) is satisfied in all cases studied 
by the author is in itself an interesting finding. That the food consump- 
tion behavior of fishes can be characterized quantitatively by an index & 
is also important. However, the author’s hypothesis is of a purely formal 
nature and does not tell us anything about the biophysical nature of the 
quantity & It would be of great interest to postulate some definite hypo- 
thetical biophysical mechanism of hunting and food consumption and to 
derive equation (2) from such a hypothetical mechanism, The mechanism 
may involve such factors as speed of locomotion of the fish, distance at 
which the fish recognizes a food particle, minimal physiologically pos- 
sible time between two successive food ingestions, etc. If equation (2) 
can be derived from such a hypothetical mechanism, then & would be ex- 
pressed in terms of the above-mentioned other quantities and could be 
calculated independently, This would considerably enhance the value 
of the agreement between expression (2) and the direct experimental data, 
An attempt in that direction, stimulated by the reading of Ivlev’s book, 
has been recently made by the reviewer, 

The author proceeds to more complicated cases and studies the effect 
of non-uniformity of distribution of food on the rate of consumption, As 
a measure of non-uniformity he introduces the quantity 


Arty 
hl": Some (3) 


where is the deviation in a given region from the average concentra- 
tion and n the number of regions. In his experiments food distributions 
are created artificially which are characterized by different values of ¢, 
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Again by very elementary and purely formal considerations the author 
atrives at a generalization of equation (2): 


r= R[1— 6 (SPted)], (4) 


where & is another constant, This equation gives the relation between 
r, the average concentration p of food, and the ‘‘index of aggregation’? 
¢. Expression (4) is found to be satisfied not only in laboratory experi- 
ments, but also in natural surroundings, namely, in the delta of the 
Volga River. The following table gives an example of the agreement in 
the latter case, 


calculated observed 
0.65 54.3 56.5 
2.33 93.6 91.5 
3.38 106.0 102.5 
1.25 61.2 56.8 
1.57 69.9 73.6 


The same remark may be made about those results as was made about 
the previously discussed ones, As a preliminary first step the purely 
formal approach is undoubtedly valuable. But further development of the 
theory will have to be made on the basis of a biophysical picture, 

Chapter IV deals with the selectivity of the food consumption, Given 
different types of available food, the fish prefers some types to others, 
The actual rates of consumption of different foods are determined both 
by the preferences and by the ‘‘accessibility’’ of the different foods, 
An appropriate index of selectivity is introduced which takes into ac- 
count the fact just mentioned. Then the author proceeds to determine 
the dependence of selectivity on various factors, notably on the state of 
hunger of the fish, These quantitative data, summarized in tables and 
curves, indicate among other things that the more hungry a fish, the 
lesser its food selectivity. To put it in common parlance, a well-fed 
fish is more choosy. This may appear to some to be a rather trivial and 
obvious conclusion which applies to humans just as well as to fishes, 
But how many times in the history of science have scientists been mis- 
led by taking for granted some ‘‘obvious’’ things! Moreover the important 
thing is that Ivlev’s data not only show that a well-fed fish is choosier, 
but they also show how much choosier it is. 

Chapter V deals with different aspects of competition for food between 
individuals of the same species, as well as between different species. 
Various aspects of the problems in this field are studied, with a defi- 
nite view to the work of Lotka, Volterra, and Gause, and with sugges- 
tions for possible generalizations. 

In the same spirit Chapter VI treats the ecology of hunger. 

This reviewer is not qualified to pass an opinion on the purely experi- 
mental side of Professor Ivlev’s work, But quite regardless of what such 
an evaluation would be, the importance of the book as contributing to- 
wards a closer cooperation of theory and experiment and as suggesting 
numerous theoretical problems will remain, 
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The reviewer is informed that two European scientists are now trans- 
lating this book into English for their personal use. It is to be hoped 
that a farsighted American publisher will use this opportunity to make it 
available to the American readers. 
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